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1.  BACKGROUND  AND  SUMMARY  OF  PHASE  I  ACCOMPLISHMENTS 

This  effort  will  incorporate  real-time  simulation  of  weapons  effects  into  a  DIS/HLA 
compatible  urban  warfare  virtual  environment.  The  Phase  I  effort  identified  key 
problems  and  relevant  approaches  for  Phase  II.  Moreover,  fundamental  modeling  and 
simulation  issues  which  extend  beyond  the  current  application  were  addressed;  i.e., 
how  to  embed  modeling  of  physical  processes  (which  is  necessary  if  the  objects  in  a 
simulation  are  to  interact  in  a  non-scripted  manner)  in  real-time  visual  simulations. 

Incorporation  of  realistic,  physically  based,  weapon  effects  into  virtual  environments 
that  can  be  used  to  train  soldiers  and  law  enforcement  personnel  represents  a  holy 
grail  for  the  military  modeling  and  simulation  community.  Currently,  visual 
representations  of  weapon  effects  are  superimposed  on  visual  simulations  almost  as  an 
“after  thought.”  As  a  consequence,  physical  processes  associated  with  employing  a 
weapon  are  not  incorporated  into  the  simulation.  For  example,  compelling  visual 
representations  of  the  muzzle  flash  and  plume  effects  from  howitzers  might  be  included 
in  a  war  gaming  virtual  environment  but  the  effects  of  cratering,  debris,  and  ejecta  (the 
purpose  of  firing  a  howitzer)  is  not  evident  on  the  characters  or  military  assets  in  the 
simulation.^  So,  while  photorealistic,  visual  animation  of  hypothetical  weapon  effects 
might  be  represented  or  “scripted”  within  a  simulation,  the  effect  on  other  objects  or 
non-scripted  events  becomes  problematic. 

At  least  three  obstacles  thwart  incorporation  of  physically  based  weapon  effects  on 
building  structural  integrity  and  personnel  in  real-time  simulations.  First,  most 
visualization  tools  merely  represent  visible  geometric  surfaces  and  not  mechanical  or 
structural  properties  necessary  to  analyze  structural  integrity  or  weapon  interaction. 
Second,  whereas  visualization  tools  effectively  represent  rigid-body  displacement 
(translation  and  rotation)  these  tools  are  limited  in  representing  non-rigid  body 
deformation  (e.g.,  crushing,  fragmentation,  and  bending).  Finally,  the  resolution 
entailed  in  describing  weapon  effects  and  their  interaction  with  objects  in  the  simulation 
is  inconsistent  with  real-time  performance  of  a  high  fidelity  visual  model. 

The  need  to  realistically  simulate  weapon  effects  with  sufficient  generality  and  degrees- 
of-freedom  in  a  training  environment  has  never  been  more  important  than  in  an  urban 
warfare  setting  where  the  rules-of-engagement,  limits  on  collateral  damage  to  non- 
combatants  and  the  urban  infrastructure,  and  weapon  performance  and  effects  can 
vary  widely. 

Collateral  effects  of  weapons,  e.g.,  fragmentation  and  backblast,  can  have  disastrous 
and  unexpected  consequences  in  these  environments  where  target  engagement 


'  Example  provided  by  D.  Reiss,  U.S.  Army  infantry  School,  personal  communication,  23  April  1996. 


M67004-96-C-0013  (2AB) 
MRC-COM-R-96<0461 


Mission  Research  Corporation 


ranges  can  be  very  close.  Urban  terrain  also  provides  tremendous  variation  in  terms  of 
weapon  effects.  For  example,  the  same  weapon  employed  against  two  different  wall 
constructions  can  have  completely  different  effects -e.g.,  ricochet  versus  penetration. 
A  5.56-mm  rifle  bullet  striking  plate  glass  windows  can  fragment,  ricochet,  or  be 
deflected  during  perforation  as  a  function  of  striking  obliquity.^  Engagement  ranges  and 
times  are  short  resulting  in  a  fast  and  confusing  pace  of  battle  where  the  risk  of  friendly 
fire,  ricochet,  fratricide,  and  injury  to  non-combatants  is  high. 

Addressing  these  issues  entails  a  high  degree  of  situational  awareness  and  mission 
planning.  This  can  only  be  realized  by  mission  rehearsal  for  all  the  unexpected 
contingencies  that  an  urban  warfare  environment  can  provide.  Further,  these 
contingencies,  particularly  as  related  to  weapon  effects,  cannot  be  safely  employed  in  a 
MOUT  training  course,  rather,  the  only  option  to  safely  demonstrate  collateral  weapon 
effects  is  in  a  virtual  training  environment.  Table  1  includes  additional  considerations 
relative  to  combat  in  built-up  areas  and  employment  of  weapons. 

A  Phase  I  SBIR  effort  that  addresses  these  issues  was  awarded  to  Mission  Research 
Corporation  (MRC)  in  February  1996  and  concluded  August  1996.  Mr.  Admiral  Piper 
was  the  U.S.  Army  Simulation,  Training,  and  Instrumentation  Command’s  (STRICOM’s) 
Contracting  Officer’s  Representative.  The  Phase  I  effort  included  eight  tasks  with  a  six 
month  period  of  performance  for  the  technical  effort. 

In  Task  1  a  scenario  was  identified  that  consisted  of  clearing  the  Town  Hall  at  the  Fort 
Benning  McKenna  MOUT  site.  We  investigated  to  what  extent  the  McKenna  Town  Hall 
represented  a  typical  structure  and  decided  to  provide  user  prescribed  options  to 
change  interior  and  exterior  construction  details.  We  also  developed  a  list  of  weapons 
and  ammunition  that  might  be  used  in  a  “precision”  MOUT  building  clearing  operation. 

In  Task  2,  we  conducted  ballistic  experiments  to  understand  the  dynamics  of  the 
interaction  of  various  weapons  with  typical  interior  building  construction.  We  also  used 
the  data  to  provide  texture  maps  for  the  Task  6  visualization  effort  and  to  develop 
analytical  models  in  the  Task  4  Structural  Integrity  effort. 

In  Task  3  we  developed  a  sub-real-time  code  which  describes  the  external  ballistics  of 
various  small  arms  of  interest.  This  code  was  favorably  correlated  with  experimental 
data  and  more  detailed  analytical  approaches.  We  also  developed  an  approach  for 
describing  the  free-field  environment  of  fragmenting  munitions  and  explosives.  For 
fragmenting  munitions  we  use  a  data  base  of  arena  test  data.  For  explosives  we  use 
modified  Sach’s  scaling  of  a  standard  blast  wave  to  obtain  the  initial  pulse  incident  on  a 
target. 


^  S.C.  Robertson,  “Rifle  Ammunition  Performance  through  Barriers,”  INTERNATIONAL  WOUND 
BALLISTICS  ASSOCIATION  REVIEW,  2(4),  1996. 
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Table  1.  Important  Characteristics  of  Combat  in  Urban  Areas^ 


a.  Hard,  smooth,  flat  surfaces  are  characteristic  of  urban  targets.  Rarely  do  rounds  impact  perpendicular 
to  these  flat  surfaces  but  at  some  angle  of  obliquity.  This  reduces  the  effect  of  a  round  and  increases  the 
threat  of  ricochets.  The  tendency  of  rounds  to  strike  glancing  blows  against  hard  surfaces  means  that  up 
to  25  percent  of  impact-fuzed  explosive  rounds  may  not  detonate  when  fired  onto  rubbled  areas. 

b.  Engagement  ranges  are  close.  Studies  and  historical  analyses  have  shown  that  only  5  percent  of  all 
targets  are  more  than  100  meters  away.  About  90  percent  of  all  targets  are  located  50  meters  or  less 
from  the  identifying  soldier.  Few  personnel  targets  will  be  visible  beyond  50  meters  and  usually  occur  at 
35  meters  or  less.  Minimum  arming  ranges  and  troop  safety  from  backblast  or  fragmentation  effects 
must  be  considered. 

c.  Engagement  times  are  short.  Enemy  personnel  present  only  fleeting  targets.  Enemy-held  buildings  or 
structures  are  normally  covered  by  fire  and  often  cannot  be  engaged  with  deliberate,  well-aimed  shots. 

d.  Depression  and  elevation  limits  for  some  weapons  create  dead  space.  Tall  buildings  form  deep 
canyons  that  are  often  safe  from  indirect  fires.  Some  weapons  can  fire  rounds  to  ricochet  behind  cover 
and  inflict  casualties.  Target  engagement  from  oblique  angles,  both  horizontal  and  vertical,  demands 
superior  marksmanship  skills. 

e.  Smoke  from  burning  buildings,  dust  from  explosions,  shadows  from  tall  buildings,  and  the  lack  of  light 
penetrating  inner  rooms  can  combine  to  reduce  visibility  and  to  increase  a  sense  of  isolation.  Added  to 
this  is  the  masking  of  fires  caused  by  rubble  and  man-made  structures.  Targets,  even  those  at  close 
range,  tend  to  be  indistinct. 

f.  Urban  fighting  often  becomes  confused  melees  with  several  small  units  attacking  on  converging  axes. 
The  risks  from  friendly  fires,  ricochets,  and  fratricide  must  be  considered  during  the  planning  phase  of 
operations  and  control  measures  continually  adjusted  to  lower  these  risks.  Soldiers  and  leaders  must 
maintain  a  sense  of  situational  awareness  and  clearly  mark  their  progress  in  accordance  with  unit 
standard  operating  procedures  to  avoid  fratricide. 

g.  Both  the  firer  and  target  may  be  inside  or  outside  buildings,  or  they  may  both  be  inside  the  same  or 
separate  buildings.  The  enclosed  nature  of  combat  in  built-up  areas  means  that  the  weapon’s  effect, 
such  as  muzzle  blast  and  backblast,  must  be  considered  as  v^ll  as  the  round's  impact  on  the  target. 

h.  Usually  the  man-made  structure  must  be  attacked  before  enemy  personnel  inside  are  attacked. 
Therefore,  weapons  and  demolition  can  be  chosen  for  employment  based  on  their  effects  against 
masonry  and  concrete  rather  than  against  enemy  personnel. 

i.  Modem  engineering  and  design  improvements  mean  that  most  large  buildings  constructed  since  World 
War  II  are  resilient  to  the  blast  effects  of  bomb  and  artillery  attack.  Even  though  modem  buildings  may 
bum  easily,  they  often  retain  their  structural  integrity  and  remain  standing.  Once  high-rise  buildings  bum 
out,  they  are  still  useful  to  the  military  and  are  almost  impossible  to  damage  further.  A  large  structure  can 
take  24  to  48  hours  to  bum  out  and  get  cool  enough  for  soldiers  to  enter. 

].  The  most  common  worldwide  building  type  is  the  12  -  to  24  -  inch  brick  building.  Table  5,  page  13,  lists 
the  frequency  of  building  types  world-wide. 


^  Source:  "An  Infantryman’s  Guide  to  Combat  In  Built  Up  Areas,"  FM  90-10-1 ,  Department  of  the  Army, 
12  Mav  1993 
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In  Task  4,  a  terminal  ballistics  code  was  developed  and  programmed  for  small  arms  of 
interest.  Code  predictions  were  correlated  with  Phase  I  ballistic  experiments  performed 
in  Task  2.  Parametric  solutions  were  programmed  for  other  terminal  ballistics 
phenomenology  (e.g.,  ricochet)  but  await  the  Phase  II  effort  for  testing.  A  stochastic 
model  was  developed  to  describe  the  distribution  of  glass  fragments  from  an  incident 
blast  wave.  Structural  integrity  for  explosive  environments  will  be  obtained  by  using  a 
commercially  marketed  hydrocode,  AUTODYN™  and  handbook  solutions  using  the 
Task  3  environments  as  a  forcing  function. 

In  Task  5,  a  methodology  was  developed  using  the  MRC  John  0.  digital  phantom.  John 
0.  will  be  interfaced  with  the  Phase  II  virtual  environment  to  describe  the  effect  of 
penetrating  wounds  from  the  projectile  environments  described  in  Tasks  3  and  4. 
Tolerance  to  overpressure/blast  effects  will  be  described  by  threshold  values  of  peak 
pressure  and  impulse  mapped  on  to  John  0. 

In  Task  6,  a  digital  “fly-through”  model  of  the  McKenna  Town  Hall  was  developed. 
Several  approaches  and  codes  were  used  to  develop  this  model;  including  Multigen™, 
Division™,  Coryphaeus™,  and  developing  entirely  new  visualization  tools.  This  was 
done  to  compare  suitability  of  these  approaches  for  visualizing  weapon  effects. 
Ultimately,  we  converged  on  a  hybrid  approach,  incorporating  “wrap-around”  code  to  an 
existing  commercial  visualization  tool.  Static  visual  representations  of  bullet  holes  were 
also  developed  to  understand  the  “choke”  points  in  real-time  visualization  of  these 
effects.  Several  approaches  were  also  developed  to  model  debris  and  ejecta. 

The  Phase  I  effort  identified  the  key  problems  and  relevant  approaches  for  Phase  II. 
Moreover,  it  attacked  fundamental  modeling  and  simulation  issues  which  extend  well 
beyond  the  current  application;  i.e.,  how  to  embed  modeling  of  physical  processes 
(which  is  necessary  if  the  objects  in  a  simulation  are  to  interact  in  a  non-scripted 
manner)  in  a  real-time  visual  simulation. 

Phase  I  Highlights 

An  important  theme  embodied  in  this  program  was  to  not  only  represent  physical 
processes  engendered  by  weapon  effects  where  everything  goes  according  to  plan  and 
deterministic  models  are  sufficient  to  describe  outcomes  but  also  to  incorporate 
analytical  models  in  sufficient  generality,  with  sufficient  degrees-of -freedom,  and 
sufficient  detail  where  things  that  can  and  do  go  wrong  can  also  be  modeled  and 
represented  in  the  simulation.  For  example,  in  Phase  I,  a  small  arms  model  was 
programmed  which  includes  deterministic  and  stochastic  modules. 

The  deterministic  small  arms  model  developed  for  the  first  phase  of  this  project 
describes  the  exfenorand  terminal  ballistics  of  rifle,  hand  gun,  and  shotgun  ammunition 
striking  various  types  of  interior  building  constructions. 
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For  rifle  and  hand  gun  ammunition,  the  non-linear  trajectory  through  air,  including  cross 
wind  effects  (which  can  be  significant  in  urban  environments  where  wind  is  channeled), 
air  drag,  and  gravity;  as  well  as  striking  velocity,  obliquity,  and  hit  location  on  the  target, 
are  described.  Projectile  striking  conditions  are  then  used  in  conjunction  with  a  library 
of  projectiles  (single  slug,  and  two  buckshot  loads  for  a  12  gauge  shot  gun,  7.62  and 
5.56  FMJ  rifle  and  9-mm  parabellum  handgun  ammunition)  and  a  library  of  interior 
building  constructions  (window,  wall,  ceiling,  floor,  and  interior  door  cross  sections)  to 
determine  subsequent  trajectory  after  interaction  with  these  targets  (i.e.,  ricochet, 
penetration,  or  perforation).  For  shot  gun  loads  with  buckshot,  in  addition  to  the 
variables  described  above,  shot  dispersion  as  a  function  of  range  is  calculated. 

The  stochastic  portion  of  the  model  which  can  be  toggled  “on/off”  with  various  options 
uses  the  deterministic  tsrget  hit  point  calculation  to  describe  the  center  of  a  random 
distribution  of  target  striking  locations.  The  distribution  about  this  center  is  described  by 
a  Circular  Error  Probability  (CEP).  The  CEP  is  obtained  from  experimental  data  bases 
which  tabulate  delivery  errors  in  terms  of  shot  dispersion  and  offset  from  an  aimpoint. 
These  data  bases,  which  are  specific  to  each  weapon  and  projectile,  have  been 
collected  by  the  Army  to  determine  the  effect  of  target  exposure  time,  training,  stress, 
and  firing  conditions  (first  shot,  single  shot,  burst  or  automatic  mode)  on  gunner 
accuracy  (see  Table  2  below  for  an  example  relative  to  target  exposure  time  using  a 
bipod  mount). 

The  external  ballistics,  deterministic  portion  of  the  model,  has  been  correlated  with  data 
supplied  by  the  Army  and  ammunition  manufactures  (see  Figure  1  and  Figure  2,  page 
10,  for  correlation  with  9-mm  M9  Pistol  Data  and  5.56-mm  M16A1  Rifle,  respectively; 
data  supplied  in  FM  101-60-27-1). 


f 


Table  2.  Delivery  Errors  for  5.56-mm  Rifle,  M16AV 


# 


# 
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■•tic  Pist«i 


Range  (ft) 


Figure  1.  Projectile  Striking  Velocity  as  a  Function  of 
Target  Range  for  the  M9,  9-mm,  Semiautomatic  Pistol 


5.5B-mn  RHIe,  M1S9A1 


Figure  2.  Projectile  Striking  Velocity  as  a  Function  of 
Target  Range  for  M16A1  Rifle  using  5.56-mm,  M193, 
Ball  Ammunition 


The  terminal  ballistics 
portion  of  the  deterministic 
model  reflects  ballistic  data 
generated  during  Phase  I 
of  this  effort.  The 
projectiles,  target  cross 
sections,  and  results  are 
shown  in  Table  3,  Figure  3, 
and  Table  4,  respectively. 

Although  no  ricochet  data 
was  available  during  the 
first  phase  of  this  effort,  a 
“first  principles”  model 
originally  developed  to 
describe  projectile  ricochet 
from  a  human  femur  in  soft 
tissue  (see  Figure  4  and 
Figure  6)  as  part  of  a 
current  DARPA/MRC  effort 
was  modified  to  describe 
ricochet  from  building 
materials.  A  parametric 
form  of  this  model  has  been 
incorporated  in  the  Phase  I 
software.  The  ricochet 
model  will  be  correlated 
with  experimental  data 
developed  during  Phase  II. 

A  comprehensive 
description  of  a  MOUT  task 
involving  common  subtasks 
of  other  MOUT  tasks;  i.e., 
clearing  a  building,  was 
developed.  In  developing 
this  task,  MRC  consulted 
with  the  U.S.  Army  Infantry 


School  (USAIS)  at  Fort  Banning,  the  FBI  Academy  in  Quantico,  Virginia,  Training 
officers  at  the  Los  Angeles  Police  Department,  and  various  Government  Training 
manuals  relative  to  Military  Operations  in  Urban  Terrain  (MOUT).  A  building  clearing 
operation  conducted  on  the  three-story  Town  Hall  building  at  the  Fort  Banning, 
McKenna  MOUT  site  was  selected  as  the  example  scenario.  The  Town  Hall  building 
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Table  3.  Ammunition  used  during  STRICOM  Ballistic  Tests  against  Targets 
Representative  of  Residential  Construction 


Gun 

Round 

Total  Mass 
(grams) 

Diameter 

(inches) 

Projectiles 
in  Salvo 

Muzzle 

Velocity 

(fps) 

Rifled  Slug 

28.4 

0.73 

1 

1600 

#4  Cu-plated 
Buckshot 

50.5 

0.24 

41 

1210 

■ISBi 

#00  Lead 
Buckshot 

29.6 

0.33 

9 

1325 

FMJ 

Parabellum 

7.5 

0.36 

1 

1155 

7.62mm 

Rifle 

FMJ 

8.1 

0.31 

1 

2100 

has  a  good  arrangement  of  rooms,  interior  walls,  stairways,  windows,  and  doors  in 
which  to  demonstrate  many  of  the  problems,  direct  and  collateral  weapon  effects, 
associated  with  using  various  weapons  in  a  MOUT  environment. 


Table  4.  Results  from  STRICOM  Phase  I  Ballistic  Experiments 


SHOT# 

TARGET 

TYPE 

ROUND 

MUZZLE 

VELOCITY 

RESIDUAL 

VELOCITY 

VELOCITY 

LOSS 

01 

INT.  WALL 

SLUG 

1600 

1500 

6% 

02 

INT.  WALL 

#00-B 

1325 

1150 

13% 

03 

INT.  WALL 

#4-B 

1210 

1090 

10% 

04 

INT.  WALL 

9mm 

1155 

1100 

5% 

05 

INT.  WALL 

7.62mm 

2100 

2080 

0% 

06 

INT.  DOOR 

#4-B 

1210 

1050 

13% 

07 

INT.  DOOR 

#00-B 

1325 

1080 

18% 

08 

INT.  DOOR 

SLUG 

1600 

1470 

8% 

09 

FLOOR 

#00-B 

1325 

980 

26% 

10 

FLOOR 

SLUG 

1600 

1480 

8% 

11 

FLOOR 

#4-B 

1210 

950 

21% 

12 

FLOOR 

9mm 

1155 

1050 

9% 

13 

CEILING 

9mm 

1155 

1090 

6% 

14 

CEILING 

#4-B 

1210 

700 

42% 

15 

CEILING 

#00-B 

1325 

650 

51% 

16 

CEILING 

9mm 

1155 

1070 

7% 

17 

CEILING 

SLUG 

1600 

1470 

8% 

18 

1/8"  STD 
GLASS 

#00-B 

1325 

N/A 

N/A 

19 

1/4"  PATE 
GASS 

#00-B 

1325 

N/A 

N/A 
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Figure  3.  Phase  I  Ballistic  Target  Cross  Sections 


Figure  4.  Frame  from  Animation 
of  Glancing  Strike  by  AK47 
Projectile  on  Idealized  Human 
Femur  Imbedded  in  Gelatin 
Matrix  (DAMD17-94-C-4099) 


Figure  5.  Elapsed  Time  in 
Microseconds  of  Glancing  Strike  by 
AK47  Projectile  on  Human  Femur 
(DAMD17-94-C-4099) 
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Table  5.  Types  of  Buildings  and  Frequency  of  Occurrence 
World-Wide 


TYPE  OF  BUILDING 

FREQUENCY  OF 
OCCURRENCE 

(Percentage) 

30-inch  Stone 

1 

8- to  10-inch  Reinforced  concrete 

6.9 

12-  to  24-inch  Brick 

63 

6-inch  wood 

16 

14-inch  steel  and  concrete  (heavy  clad) 

2 

7-inch  steel  and  concrete  (light  clad) 

12 

The  initial  set  of  Phase  I  domestic  weapons  considered  was; 

•  M16A1,  .22  caliber  rifle  firing  5.56-mm,  M193  Ball  rounds; 

•  M16A2,  .22  caliber  rifle  firing  5.56-mm,  M855  Ball  rounds; 

•  M249,  Squad  Automatic  Weapon  (SAW),  22  caliber  machine  gun  firing  5.56-mm, 
M855  Ball  rounds; 

•  M60,  .30  caliber  machine  gun  firing  7.62-mm,  M80  ball  rounds 

•  9-mm  Semiautomatic  M9  pistol  firing  M822  ball  rounds 

•  12-gauge  “sawed-off’  shotgun  firing  single  slug,  00  and  number  4  buckshot  (not 
standard  military  issue);  and, 

•  40-mm  M203  Grenade  Launcher  firing  M406  HE  and  M433  HEDP  rounds. 

Three  soldiers  composed  a  fire  team,  and  at  least  two  fire  teams  were  involved  in  a 
clearing  operation. 

It  was  observed  that  construction  details  of  typical  buildings,  as  well  as  corresponding 
weapon  effects,  can  vary  significantly  from  the  cinder  block  and  poured  concrete 
construction  of  the  McKenna  Town  Hall  building  (see  Table  5).  For  this  reason 
alternative,  but  more  typical,  wall,  ceiling,  floor,  glazing,  and  door  constructions  were 
considered  in  Phase  I.  In  Phase  II,  a  library  of  user  prescribed  construction  details  will 
be  included  so  that  construction  details  of  building  components  can  be  varied  by  the 
user  during  the  simulation.  The  software  will  prompt  the  user  for  these  details  which 
should  be  the  same  considerations  that  occur  to  the  user  before  employing  weapons  in 
a  real  MOUT  mission.  Weapon  effects  algorithms  will  describe  differences  in  effects 
due  to  variation  of  construction  methods  and  materials. 

During  Phase  I  Ballistic  tests  employing  weapons  and  ammunition  above  on  typical 
interior  wall,  ceiling,  floor,  door,  and  glazing  materials  were  conducted  at  a  range  of 
eight  feet.  Exit  velocity  of  the  striking  projectile  through  the  rear  surface  of  the  target 
was  determined,  debris  fields  were  video  taped,  ejecta  was  captured  in  witness  plates, 
and  the  target  entrance  and  exit  holes  were  photographed.  These  tests  were  conducted 
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Table  6.  Real  Time  versus  Execution  time  for  Software 
Algorithms  in  Small  Arms  Model 


Muzzle 

Range 

Elapsed  time 

Execution 

Velocity 

(ft) 

to  strike 

time  on  66 

(fps) 

target  after 

Mhz  80486  PC 

weapon  fired 

(sec) 

(sec) 

1000 

10 

0.01 

1000 

100 

0.02 

2000 

10 

0.0052 

0.01 

2000 

100 

0.052 

0.01 

to  capture  phenomenology,  yield  preliminary  data  on  projec^le  velocity  loss  after  target 
penetration,  and  provide  texture  maps  for  the  Task  6  Visualization  and  Scene 
Animation  task. 

Software,  which  executes  in  sub-real-time  (see  Table  6)  was  also  developed  to 
describe  the  exfer/orand  terminal  ballistics  (including  the  effect  of  cross-wind,  ricochet, 
and  penetration)  for  handgun,  rifle,  and  shotgun  fired  projectiles  of  interest  (shot 
dispersion  as  a  function  of  range  is  also  described  for  shotgun  rounds).  Single  shot, 
burst,  and  automatic  modes  were  included  in  the  model  as  well  as  the  effect  on 
aimpoint  accuracy  (in  terms  of  circular  error  probability  as  a  function  of  range)  of. 


•  target  exposure  time; 

•  soldier  stress; 

•  training;  and, 

•  typical  shot  dispersion  relative  to  an  aimpoint. 

A  stochastic  model  was  also  developed  to  describe  distribution  of  glass  shards  (in 
terms  of  ejecta  mass  and  velocity)  from  small  arms  and  blast  waves.  An  approach 
which  will  be  implemented  in  Phase  II  (described  in  the  Phase  II  proposal)  was  also 
developed  to  describe  the  initial  shock,  hole  size,  and  time  and  spatial  evolution  of 
overpressure  from  the  blast  and  fragmentation  from  grenades  and  other  larger 
ordnance  (such  as  demolition  satchels). 

The  effect  of  penetrating  wounds  on  personnel  is  accomplished  using  the  “John  O.” 
phantom  developed  by  MRC  in  a  U.S.  Army  Natick  Phase  II  SBIR  (DAAK60-C-92-008) 
-  see  Figure  6,  Figure  7,  Figure  8,  and  Figure  9.  This  digital  representation  of  a  soldier 
divides  the  human  anatomy  into  80,000  volume  elements  and  250  tissue  types.  The 
anatomy  can  be  articulated  and  scaled  to  vary  anthropometry.  Currently,  body  armor 
and  ballistic  protective  goggles  can  be  placed  on  John  0.,  dynamic  exposure 
probabilities  of  user  prescribed  regions  on  the  human  body  during  task  performance 


Figure  6.  Autopsy  Mode  of 
"John  O.”  (MRCMAN)  Phantom 


Figure  7.  Projection 
Mode  of  John  O 


u 


can  be  obtained,  and  the  effect  of  penetrating  wounds  from  virtual  small  arms  can  be 
demonstrated  on  John  0.  As  of  September  1996,  John  0.  has  resulted  in  almost  $2M 
of  Phase  III  contract  awards. 

A  digital  representation  of  the  McKenna  Site  Town  Hall,  that  can  be  explored  in  real¬ 
time  using  a  Silicon  Graphics  workstation  was  also  developed.  This  was  done  to 
identify  constraints  associated  with  incorporating  real-time  representations  of  weapon 
effects  in  the  Phase  ll  virtual  environment.  A  narrated  video  tape  of  a  user  exploring  the 
MRC  Digital  Town  Hall  in  real-time  was  sent  under  separate  cover  to  STRICOM 
(Admiral  Piper). 
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2.  SMALL  ARMS  EXTERIOR  AND  TERMINAL  BALLISITCS  MODEL 


This  section  discusses  the  analytical  and  numerical  background  of  the  MRC  developed 
computer  code,  MRSAW  (Mission  Research  Small  Arm  Weapons).  A  listing  of  the 
MRSAW  source  code  is  presented  in  Section  5  of  this  report.  This  code  was  developed 
during  the  first  phase  of  this  effort  and  is  capable  of  performing  the  analyses  listed 
below  in  sub-real-time.  A  Nomenclature  describing  all  variables  and  parameters  used 
in  Section  2  is  presented  in  Table  7,  page  17. 


I.  EXTERIOR  BALLISITCS.  Motion  of  bullets  in  air  prior  to  hitting  the 

TARGET  UNDER  PRESCRIBED  INITIAL  CONDITIONS,  I.E.,  MUZZLE  VELOCITY,  ANGLE  MADE 
BY  THE  INITIAL  BULLET  PATH  WITH  THE  HORIZONTAL,  AND/  OR  TARGET  LINE-OF-SIGHT 

II.  TERMINAL  BALLISITCS 

In  the  case  of  target  penetration.  MRSAW  calculates 

•  Projectile  striking  velocity  on  the  target 

•  Projectile  exit  velocity  from  the  target 

•  T OTAL  energy  deposited  ON  THE  TARGET  (RELATED  TO  CREATION  OF  DAMAMGE) 

•  Damage  area  of  the  target  in  the  case  of  buckshot  or  slugger  rounds 

FIRED  FROM  A  SHOTGUN 

In  the  case  of  no  penetration.  MRSAW  calculates 

•  Bullet  ricochet 

•  PROJECTILE  MOTION  AFTER  REBOUND 


Target  locations  and  material  properties  need  to  be  specified  by  the  user  in  order  to 
implement  calculations  in  items  (II)  and  (III), 

The  first  phase  of  projectile  motion,  referred  to  as  projectile  exterior  ballistics,  describes 
the  projectile  trajectory,  prior  to  impact,  subject  to  gravity,  aerodynamic  drag,  and  cross 
wind  effects.  Projectiles  fired  from  rifles  and  other  small  arms  as  well  as  projectiles 
launched  from  fragmenting  munitions  (grenades,  artillery  shells,  and  bombs)  are 
considered.  If  aerodynamic  drag  can  be  specified  exactly,  the  problem  of  determining 
the  projectile  trajectory  is  deterministic  even  though  the  resulting  equations  of  motions 
may  be  nonlinear.  Under  these  circumstances,  an  exact  closed  form  solution  for  the 
projectile  motion  is  not  possible. 
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Table  7.  Section  2  Nomenclature 


# 


# 


# 


# 


Symbol 

Meaning 

Where  defined 

A 

Velocity  modified  drag  coefficient 

Page  23,  equation  (6a) 

a 

acceleration 

Page  20,  equation  (1) 

C(SJ 

Ballistic  coefficient 

Page  20 

D 

Lateral  projectile  dimension 

Page  31 

A 

Projectile  diameter  in  inches 

Page  21 ,  equation  (3) 

d 

Target  thickness 

Page  29 

Ear 

Minimum  energy  needed  for  complete  target 
penetration  for  a  given  projectile-target  pair 

Page  30 

e 

Coefficient  of  restitution  between  a  given 
projectile-target  pair 

Page  33 

Fo 

Defined  constant  (=— ) 

A 

Page  24 

Fr 

Constant  force  of  retardation  during  target 
penetration 

Page  29 

9 

Acceleration  due  to  gravity 

Page  23 

L 

Longitudinal  projectile  dimension 

Page  31 

N 

Total  number  of  pellets  in  a  shot 

Page  34 

m 

Projectile  mass 

Page  29,  equation  (12) 

n 

=  n{v)  from  Ingall’s  Table  8 

Page  21 

n(v) 

Power  index;  Mayevski’s  parameter 

Page  20,  equation  (2) 

n(R,r) 

Number  of  pellets  within  radius  rat  range  R 

Page  33 

nl,n2 

Experimental  data  points  for  n{R,r) 

Page  34 

R 

Range  in  feet 

Table  10,  Page  25 

R1,R2 

Experimental  points  for  range  R 

Page  34 

Ff 

Reference  range  for  experimental  data 

Page  34,  equation  (15) 

r 

Radial  coordinate  on  the  target  plane 

r 

max 

Spray  radius  for  damage  on  the  target 

Page  39 

A 

Shape  class  parameter 

Page  21 

V 

Projectile  speed 

Page  19 

^0 

Muzzle  velocity 

Page  19 

^in 

Incident  projectile  velocity  on  target 

Page  29,  equation  (13) 

^ex 

Projectile  exit  velocity  after  target  penetration 

Page  29,  equation  (13) 

Vcr 

Minimum  incident  velocity  needed  for  complete 
target  penetration  for  a  given  projectile-target  pair 

Page  30 

W, 

Projectile  mass  in  grains 

Page  21,  equation  (3) 
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Nomenclature-Continued 


Symbol 

Meaning 

Where  defined 

x(t) 

Horizontal  location  of  the  projectile  at  time  t 

Page  22 

X 

Projectile  speed  along  the  x-direction 

Page  23 

Vertical  location  of  the  projectile  at  time  t 

Page  22 

HBH 

Projectile  speed  along  the  y-direction 

Page  24 

z(t) 

Lateral  location  of  the  projectile  at  time  t 

Page  22 

z 

Projectile  speed  along  the  y-direction 

Page  24 

a 

overhead  dot  indicates  time  derivative  of  any  a 

Section  2 

a(v) 

Meyevski’s  parameter 

Page  20,  equation  (2) 

m 

Air  drag  factor  due  to  compression 

Page  20,  equation  (1) 

5 

Shooting  angle  with  the  horizontal 

Page  22,  Figure  10 

X 

Model  parameter  for  n{R,r) 

Page  34,  equation  (15) 

A 

Model  parameter  for  n(jR,r) 

Page  34,  equation  (15) 

9 

Spray  angle  for  pellets 

Page  34 

2.1  Rationale  Behind  Using  the  Full,  Nonlinear,  Coupled  Theory  to  Describe 
Projectile  Motion 

When  a  projectile  is  launched  in  space  either  as  a  bullet  from  a  gun  or  as  a  fragment 
from  a  blast,  its  trajectory  in  space  depends  on  the  following  geometric  and  kinematic 
properties; 

1.  Initial  speed  at  the  time  of  its  release  in  air; 

2.  Initial  direction  of  its  velocity  vector  at  the  time  of  its  release  in  air; 

3.  Air  drag  coefficient  which  is  the  constant  of  proportionality  in  the  direct  variation  of 
projectile  deceleration  with  the  square  of  the  instantaneous  projectile  speed  (This  is 
in  accordance  with  the  Newton’s  law  of  aerodynamic  drag);  and, 

4.  Gravitational  acceleration  due  to  the  weight  of  the  projectile. 

Size  and  shape  control  rotational  motion  about  the  projectile  center  of  mass  while  the 
parameters  listed  above  (items  1  through  4)  control  the  trajectory  of  the  projectile 
center  of  mass. 

For  cases  where  the  initial  velocity  of  the  projectile  is  very  high  (more  than  100  ft/sec), 
and  short  flight  paths  are  considered  (50  to  100  ft),  gravitational  effects  can  be  ignored 
so  that  the  projectile,  idealized  as  moving  in  a  straight  line,  yields  acceptable  results. 
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Strictly  speaking,  under  any  circumstances,  when  gravity  is  present  and  the  projectile  is 
not  moving  along  a  vertical  line,  the  trajectory  of  the  center  of  mass  of  a  projectile  is 

curved  so  that  its  curvature,  - ,  mass,  m ,  and  speed,  v ,  is  related  to  the  normal 

P 

component,  F„ ,  of  forces  resulting  from  the  air  drag  and  gravity  according  to  Newton’s 

second  law  of  motion,  viz.,  F„  =  m— .  According  to  this  equation,  the  issue  of  whether 

P 

the  assumption  of  a  straight  line  path  of  the  projectile  is  valid  or  not,  really  hinges  on 

F„ .  This  result  can  be  put  in  a  better  perspective  by  writing  this  equation  as  —  = 

p  mv 

where  a  is  a  characteristic  dimension  of  the  projectile.  Thus,  curvature  can  be  ignored 
F  a 

when  « 1 .  With  decreasing  projectile  velocity,  this  condition  ceases  to  be  valid 
mv 

after  some  point  in  time  (see  Figure  10). 

In  order  to  develop  a  code  which  is  capable  of  predicting  projectile  location  and  velocity 
at  any  time  after  it  is  launched  in  space,  we  can  not  in  general  depend  on  a  model 
which  relies  on  a  linear  projectile  path.  That  is,  we  want  the  ability  to  describe 
projectiles  launched  from  building  interiors  to  the  outside  as  well  as  the  ability  to 
describe  projectiles  launched  from  the  outside  into  building  interiors.  We  also  want  the 
ability  to  model  collateral  effects  which  may  involve  debris  and  fragments  ejected  over 
large  distances.  These  situations  may  involve  long  flight  paths  as  well  as  large  forces 
retarding  projectile  motion  which  violate  the  condition  above.  In  the  absence  of 
channeled  air  flows,  linear  assumptions  are  probably  acceptable  for  projectiles 
launched  within  building  interiors  at  interior  targets  since  these  involve  high  projectile 
velocities  and  short  flight  paths. 

If  we  assume  that  the  aerodynamic  drag  on  the  projectile  is  known  as  function  of  the 
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instantaneous  projectile  speed,  then  the  resulting  equations  of  motion,  in  the  presence 
of  gravity,  yield  three  (for  full  3-D  problem)  or  two  (for  a  two-dimensional  problem) 
coupled,  nonlinear  second  order  differential  equations.  These  equations  can  be 
numerically  solved  by  using  a  modified  version  of  the  Runge-Kutta  intergration 
algorithm  for  solving  nonlinear,  second  order,  differential  equations. 

For  high  projectile  velocities  and  short  flight  paths,  these  equations  reduce  to 
uncoupled,  nonlinear,  second  order  differential  equations.  These  equations  possess 
complete  analytical  solutions.  The  numerical  code  developed  under  this  program  yields 
numerically  integrated  results  from  the  coupled  theory  and  analytical  solutions  from  the 
uncoupled  theory.  For  most  applications,  execution  times  to  obtain  both  solutions  are 
comparable.  However,  if  results  from  the  uncoupled  equations  are  acceptable,  then  this 
solution  should  be  used  when  the  code  is  used  repeatedly  to  analyze  a  series  of  these 
types  of  problems  (e.g.,  grenade  fragments  and  automatic  rifle  fire  inside  a  room). 

According  to  conventional  formulations  of  aerodynamic  drag,  the  force  of  retardation  or 
drag  acting  on  a  projectile  moving  in  air  or  fluid  with  a  speed,  v ,  is  directly  proportional 
to  the  projected  cross-sectional  area  of  the  projectile  on  a  plane  normal  to  the  direction 
of  the  velocity  of  the  projectile  and  the  square  of  the  speed,  v ,  i.e.,  v^.  However, 
depending  on  the  range  of  v ,  the  constant  of  proportionality  may  depend  on  the  speed, 
V .  Using  the  above  law  of  aerodynamics,  the  deceleration  -a  of  the  projectile  can  be 
written  as: 


-a  =  /3(vy\  (1) 

where  p(y)  is  related  to  the  Meyevski’s  parameter  a{v)  and  index  n{v)  through  the 
equation; 


m= 


C(Se) 


(2) 


In  Equation  (2),  C(S^)  is  the  ballistic  coefficient  which  is  a  function  of  the  bullet  shape 
class,  .  The  concept  of  a  bullet  shape  class  that  influences  retardation  of  a  projectile 

in  air  stems  from  development  of  a  reference  projectile  with  which  other  projectiles  are 
scaled.  For  our  purposes,  the  reference  projectile  has  a  weight  of  one  pound,  a 
diameter  of  one  inch,  and  a  ballistic  coefficient,  C  =  1  For  the  reference  projectile,  the 
dependency  of  the  Meyevski’s  parameter,  a{v) ,  and  index,  n{v),  on  the  instantaneous 
speed,  V ,  of  the  projectile  is  shown  in  Ingall’s  table  (Table  8). 

Ingall’s  table  indicates  that  for  low  projectile  velocities  {v<  790 fps )  n=0  so  that  the 
retardation  formula  presented  in  Equation  (1)  reduces  to  the  classical  Newton 


k.  J.  Pejsa,  MODERN  PRACTICAL  BALLISTICS,  Second  Edition,  Kenwood  Publishing,  1989. 
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9  Table  8.  Ingall's  Table:  Retardation  coefficients  a(v)  and  n(v)  for  various 

velocities,  v,  for  the  reference  projectile. 


a(v) 

n(v) 

V  (fps) 

4.065  E-3 

0.45 

2600  <  V  <  3600 

1.248  E-3 

0.30 

1800<  v<2600 

0.01316  E-3 

0 

1370<  v<  1800 

0.957  E-7 

-1 

1230<  v<  1370 

0.6337  E-13 

-3 

970  <  V  <  1230 

0.5935  E-7 

-1 

790  <v<  970 

0.4676  E-4 

0 

0<v<  790 

Table  9.  Shape  Class  Descriptions 


Shape  Class  (Sc) 

Projectile  Shape 

1 

Slender  Spitzer  Boat  tail 

2 

Flat  Base  Spitzer 

3 

Semi-round  nose  or  semi-Spitzer 

4 

Round  nose  or  flat  nose  common  in  30-30 
cartridge 

Note:  In-between  shape  class  may  also  be  used  in  M 

RCSAW  code. 

formulation  of  aerodynamic  drag  with  a  dependency.  For  higher  velocities,  the 
power  index  of  v  as  well  as  Meyevski’s  parameter  depends  on  projectile  velocity  due  to 
increasing  air  compression  (air  density)  in  front  of  the  projectile  with  increasing 
projectile  velocity. 

The  algorithm  used  in  MRSAW  code  to  calculate  the  parameters  a(v),  n(v),  and  C(Sc) 
is  as  follows.  Ingall’s  table  is  included  in  the  database  as  a  lookup  table.  For  any  input 
muzzle  velocity,  these  parameters  are  read  from  the  table,  and  are  used  in  the  drag 
calculation.  For  the  determination  of  the  ballistic  coefficient,  C(SJ ,  for  a  given 

projectile,  we  use  a  shape  class  parameter.  Sc,  which  is  an  index  parameter  that 
describes  the  shape  of  the  projectile.  The  index  for  the  shape  parameter  Sc  is  shown  in 
Table  9. 


Using  a  well  established  empirical  relation,  the  indexed  shape  parameter  Sc  is  used  to 
calculate  the  ballistic  coefficient,  C,  for  any  given  projectile.  This  relation  is  given  by 


C(Sc)  = 


%+0.5'^ 


415(8-2Sc+Sc  ) 


(3) 
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In  equation  (3),  is  the  projectile  mass  in  grains  (1grain=0.0648  gram)  and  D,  is  the 
projectile  diameter  in  inches. 


2.2  Determination  of  the  First  Phase  of  Motion  in  Exterior  Ballistics:  Equation  of 
Motion  of  a  projectiie  moving  through  air 


The  kinematics  of  a  projectile  moving  in  air  is  shown  in  Figure  1 1 .  Assuming  that  the 
bullet  of  mass  m  is  moving  in  the  xy-plane,  Newton’s  second  law  of  motion  yields  the 
following  equations  of  motion  along  the  x-  and  y-axis  respectively. 


X  = 


-a{v)v 


C(S^) 


y  =  -9- 


c(Sc)  yj 


(4a) 

(4b) 


;  z(t) =0 

In  Equation  (4),  we  have  used  overhead  dots  to  indicate  time  derivatives. 

From  the  nature  of  the  differential  equations  in  (4),  it  is  obvious  that  they  are  nonlinear 
and  hence  can  not  be  integrated  in  closed  form  to  express  the  velocity  and  location  of 
the  projectile  as  a  function  of  time. 

2.3  Linearization  of  Equation  (4) 

Under  some  restrictive  cases,  it  is  possible  to  linearize  Equation  (4)  so  that  a  closed 
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form  solution  for  the  projectile  in  the  hodograph  (locus  of  speed  as  a  function  of 
location)  in  the  physical  plane  can  be  derived.  The  rationale  behind  such  linearization 
is  discussed  below. 

For  bullets  shot  at  high  speeds  (  >  1000  fps)  at  low  angles  (  <  5  degrees)  for  short 
ranges  (  <  100  ft.),  the  bullet  path  is  mostly  horizontal  so  that  the  velocity  along  the 
horizontal  direction  is  much  larger  than  the  velocity  in  the  vertical  direction.  Besides,  no 
significant  reduction  in  velocity  occurs  during  the  travel  time.  In  view  of  these 
considerations,  the  following  assumptions  can  be  justified. 

With  reference  to  Figure  1 1 ,  we  assume  that  the  bullet  location  at  any  time  t  is  given  by 
the  coordinates  {x{t),y(t),z{t)}  where,  in  the  absence  of  any  cross-wind  in  the  z- 
direction,  z(t)  =  0.  At  the  end  of  this  section,  we  will  discuss  the  effect  of  cross-wind 
along  the  z-direction.  The  exact  equations  of  motion  (4)  can  be  rewritten  as  [Reference. 
4]. 


x=  —  =  V— =  1  + 

dt  dx  C(Sc)  t  x'; 


-a(y)  ^.2-n(t>) 

C(S^) 


1 

v  =  {x^+y^y 


X 


2  . 


for 


{5a) 


y  =  -9- 


a{v) 


yv 


\-n 


(5b) 


Assuming  that  in  the  range  of  v  considered,  a(v)  and  n(v)  are  constants  so  that  for  the 
reference  projectile  {C=1),  Equation  (5)  becomes 


x=^  =  v—  =  -Av^-^-,  where  A  =  a(v)/C(Sc) 
dt  dx 


(6a) 


y  =  -g-Ayv'^’'  (6b) 

In  Equations  (6),  a{v)  =  constant  =  A  and  n(v)  =  constant  =  n  whose  values  can  be 
computed  from  Ingall’s  Table,  Table  8,  for  a  reference  projectile. 

Integrating  (6a),  yields; 

x  =  v=Vo(\-—Y  (7) 

r  n 
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In  Equation  (7),  we  have  used  the  initial  condition  that  x(0)  =  Vq  Fq 

Defining  a  function  F(t)  =  Ay"” ,  we  see  that  F(0)  =  and  F  =  F^-nx 
Integrating  (7)  again  with  respect  to  f,  yields: 


x=- 


n 


1-1  + 


(l-n)yot  " 


(8) 


Similarly,  integrating  (6b),  we  have 

xqFq  if  nx 


y  = 


{n-2)v^ 


1— 


1- 
V  ^’oy 


-1 


(9) 


where  we  have  used  the  initial  condition  y(0)  =  0  (small  angle  approximation). 
Integration  of  Equation  (9)  again,  yields: 


y  = 


gFo 


(n-2)Vo 


2(1 -n) 


v2— 


^  nx 


-X 


(10) 


In  MRSAW,  we  have  used  equations  (7)-(10),  to  calculate  the  velocity  and  location  of 
the  projectile  at  any  time  t  (or  any  range  x=R)  when  the  initial  muzzle  conditions  are 
given.  This  solution  has  been  referred  to  as  ‘Linear  Theory  Results’  in  MRSAW. 

2.4  Some  Drawbacks  of  the  Linearized  Solution 

From  the  initial  conditions  used  in  deriving  the  linearized  solutions  (7)-(10),  it  is  obvious 
that  the  solution  is  independent  of  the  initial  angle  of  shooting  S  with  the  horizontal 
(Figure  11).  As  this  angle  significantly  influences  the  rise  and  drop  of  the  bullet  from  the 
horizontal  line,  solutions  presented  by  equations  (7)-(10)  can  not  be  expected  to  be 
accurate  in  estimating  the  rise  and  drop  of  the  bullet  for  any  given  S.  This  can  be 
verified  later  when  the  solution  from  the  exact  equations  are  presented.  Besides,  the 
basis  of  the  linearization  process  is  dependent  on  the  accuracy  of  the  condition  that 
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Table  10.  Comparison  between  linearized  solution  and  exacy  solution  for  the 
reference  projectile  with  an  initial  velocity  of  3,000  fps 


Angle,  5 
(Degrees) 

Range,  R 
(feet) 

Exact  Soln. 

For  Impact 
Velocity  (fps) 

Linear  Soln. 

For  Impact 
Velocity  (fps) 

Exact  Soln. 

Rise  (+)  or 

Drop  (-)  in 
inches 

0 

10 

2997 

2997 

-  0.002 

0 

25 

2992 

2992 

-  0.0124 

0 

50 

2984 

2984 

-  0.054 

0 

100 

2968 

2968 

-  2.15 

0 

1000 

2687 

2687 

-  22.57 

1 

10 

2997 

2997 

2.02 

1 

25 

2992 

2992 

5.02 

1 

50 

2984 

2984 

10.41 

1 

100 

2967 

2968 

20.65 

1 

1000 

2686 

2687 

184.32 

2 

10 

2997 

2997 

4.03 

2 

25 

2992 

2992 

10.05 

2 

50 

2984 

2984 

20.86 

2 

100 

2967 

2968 

41.52 

2 

1000 

2686 

2687 

391.20 

5 

10 

2997 

2997 

10.49 

5 

25 

2992 

2992 

26.15 

5 

50 

2984 

2984 

52.20 

5 

100 

2967 

2968 

104.0 

5 

1000 

2686 

2687 

1010.76 

Linear  Soln 
Rise  (+)  or 
Drop  (-)  in 
inches 


-  0.002 

-  0.0124 

-  0.054 

-  2.15 

-  22.48 


2.01 

5.02 

10.41 

20.66 

184.08 


.03 
.05 
.88 
.52 
390.94 


9 


104.41 

1012.92 


« 1  which  becomes  a  bad  approximation  as  the  shooting  angle  S  with  horizontal 


increases  or  when  the  range  increases. 

For  applications  in  buildings,  range  is  not  really  an  issue  since  the  target  range  is 
below  100  ft.  for  most  applications,  but  the  shooting  angle  can  be  arbitrary  since  either 
the  wall,  or  the  floor  or  the  ceiling  may  be  targets. 


2.5  Exact  Solution  of  the  Projectile  Equation  of  Motion 

The  exact  solution  of  the  problem  of  describing  the  trajectory  of  the  bullet  can  be 
obtained  from  the  differential  Equation  (4)  under  the  general  initial  conditions  (Figure 
11): 


x(0)  =  Vq  cosS 
y{0)  =  V(^  sin^ 
xi0)  =  0  =  y(0) 
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Since  the  differential  equations  in  (4)  are  nonlinear,  no  closed  form  solution  can  be 
obtained  as  in  the  case  of  the  linearized  form  of  these  equations  shown  in  (5).  Solution 
has  been  obtained  from  an  extension  of  the  highly  accurate  algorithm  for  solving  one¬ 
dimensional,  nonlinear,  second  order  differential  equations  known  as  the  Runge-Kutta 
method  (RKM).  The  extended  version  of  this  algorithm  can  solve  either  two  (for  planar 
projectile  motion)  or  three  (for  fully  three-dimensional  projectile  motion)  coupled 
nonlinear  differential  equations.  The  validity  of  this  extended  version  of  RKM  has  been 
verified  by  comparing  the  solution  from  RKM  with  exact  solutions  for  coupled,  second 
order  differential  equations  where  exact  solutions  are  available.  This  version  of  the 
solution  is  referred  to  as  ‘exact  solution’  in  MRSAW. 

Some  effect  of  the  initial  angle  6  on  the  linearized  solutions  (7)-(10)  can  be 
incorporated  as  follows.  The  rise  of  the  bullet  in  the  absence  of  any  vertical 
deceleration  due  to  gravity  is  given  by  y  =  xtan(<5)  so  Equation  (10)  can  be  modified  by 
the  addition  of  this  term  to  the  right  hand  side  of  (10).  This  modification  has  been 
included  in  MRSAW. 

Some  comparison  of  the  results  predicted  by  the  exact  and  linearized  solution  of  the 
bullet  flight  path  can  be  made  at  this  point.  Let  us  consider  a  target  at  distance  varying 
between  10  and  1,000  feet.  For  the  reference  bullet,  results  for  a  muzzle  velocity  of 
3,000  fps  and  various  shooting  angles  are  presented  in  Table  10.  It  is  clear  from  Table 
10  that  results  are  in  good  agreement  for  horizontal  location  of  the  bullet  while  vertical 
location  predicted  by  the  linear  solutions  are  in  error  for  shooting  angles  above  4 
degrees. 

For  any  application  results  from  both  types  of  solutions  are  presented.  The  logic  to 
automatically  select  the  appropriate  solution  will  be  included  in  the  Phase  II  software. 
The  rest  of  the  MRSAW  development  described  below  employs  the  exact  solution  form 
of  the  equations  of  motion  shown  (Equations  4). 

2.6  Second  Phase  of  the  Exterior  Ballistics 

The  first  phase  of  the  exterior  ballistics  deals  with  the  flight  path  of  the  projectile  from 
the  weapon  to  the  target.  At  the  time  of  hitting  the  target,  MRSAW  presents  the  data 
just  before  target  impact.  These  data  include:  (1)  Projectile  striking  velocity;  and,  (2) 
Projectile  impact  angle  on  the  target. 

What  happens  after  the  projectile  hits  the  target  depends  on  the  material  and 
geometrical  properties  of  the  projectile-target  pair.  Following  are  the  possible 
scenarios. 


In  the  case  of  a  bullet,  we  may  have: 
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1 .  Bullet  incapable  of  penetrating  the  target  and  the  force  of  impact  may  cause  the 
bullet  to  ricochet  after  impact; 

2.  If  ricochet  occurs,  the  plane  of  motion  of  the  bullet  may  differ  from  the  plane  of 
motion  (for  two-dimensional  motion)  prior  to  impact.  The  plane  of  motion  after 
impact  is  the  plane  that  contains  the  bullet  line  of  impact  and  the  line  offeree  of 
impact  from  the  target.  In  MRSAW,  the  line  offeree  of  impact  between  the  projectile 
and  the  target  is  assumed  to  be  normal  to  the  target  at  the  point  of  impact.  For 
various  projectile-target  pairs  of  interest,  the  nature  of  the  line  of  impact  needs  to  be 
determined  from  experiments;  and, 

3.  If  the  bullet  penetrates,  the  exit  velocity  of  the  bullet  is  of  interest  (Terminal  Ballistic 
phase). 

In  the  case  of  buckshot  evaluated  during  Phase  I  and  launched  from  shot  guns, 
ricochet  does  not  occur  (since  the  buckshot  is  made  predominately  from  lead  and  will 
therefore  tend  to  deform).  These  projectiles  may  however: 

1 .  Rebound  after  hitting  the  target.  In  this  case,  its  flight  path  after  impact  needs  to  be 
determined;  and, 

2.  Penetrate  the  target  in  which  case  its  exit  velocity  is  of  interest. 

Since  the  solution  for  the  post  impact  conditions  depends  on  the  projectile-target  pair 
and  their  properties,  we  have  used  test  data  for  various  projectiles-targets 
combinations.  MRSAW  uses  these  test  results  to  predict  the  exit  velocity  when  the 
projectile  penetrates  the  target.  In  the  case  of  projectile  rebound,  the  force  of  impact  is 
calculated  from  an  assumed  coefficient  of  restitution  which  also  needs  to  be  determined 
from  experiments  during  Phase  II  for  various  projectile  target-pairs.  For  the 
determination  of  the  second  phase  of  exterior  ballistics,  MRSAW  has  the  following 
projectiles  and  target  properties  in  its  data  base.  This  data  base  is  derived  from  the 
MRC  experiments  during  the  first  phase  of  this  program,  and  is  presented  in  a  tabular 
from  in  Table  11. 

The  library  of  projectiles  currently  included  in  MRSAW  are  listed  below. 

1 .  Shotgun  Single  Slug 

2.  Shotgun  #00  lead  spheres 

3.  Shotgun  #04  copper  plated  lead  spheres 

4.  9-mm  parabellum  handgun  bullet 

5.  7.62-mm  and  5.56-mm  rifle  bullets 

The  library  of  targets  currently  included  in  MRSAW  are  listed  below. 

1 .  Generic  interior  wall 

2.  Generic  interior  floor 

3.  Generic  interior  door 

4.  Generic  ceiling 
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The  properties  of  these  projectiles  and  targets  are  given  in  Table  3,  and  Figure  3. 

In  the  case  of  penetration,  exit  velocity  has  been  calculated  using  the  results  in  Table 
1 1  as  follows. 

Assuming  that  the  retarding  force,  ,  during  penetration  is  a  constant,  the  amount  of 
energy  loss  due  to  the  penetration  is  equal  to  the  product  of  F^.d  where  d  is  the  depth 
of  the  target  in  the  direction  of  penetration.  This  result  can  be  shown  by  integrating  the 
one-dimensional  equation-of-motion  (assuming  that  the  penetration  process  can  be 
modeled  as  an  equivalent  one  dimensional  problem)  of  a  projectile  of  mass,  m, 

mx—  =  mx=^-F  (12) 

dx  ^ 

Integrating  Equation  (12)  from  x  =  0  to  x  =  d,  we  have 

-vJ)  =  F4  (13) 

In  Equation  (13),  the  subscripts  ‘in’  and  ‘ex’  refer  to  quantities  just  before  and  after 
penetration  of  the  target.  Since  the  left  hand  of  Equation  (13)  represents  energy  loss, 
AE  ,  due  to  penetration.  Equation  (13)  indicates  that  column  10  of  Table  1 1  is  equal  to 
the  product  F^.d  of  the  average  force  of  retardation,  F^ ,  and  the  depth  of  penetration, 
d .  Assuming  that  F^.d  is  a  property  of  the  projectile-target  pair.  Equation  (13)  can  be 
rewritten  as  Equation  (14)  to  yield  the  exit  velocity,  ,  when  the  entry  velocity,  ,  is 
given. 


V 


ex 


1 


F4 


1  2 

-mVi„ 

2 


2.7  Determination  of  Buiiet  Ricochet  and  Siug  Rebound 


(14) 


In  the  context  of  this  report,  ricochet  refers  to  motion  of  the  projectile  which  is  described 
in  terms  of  both  motion  of  the  projectile  center-of-mass  and  rotation  about  the  projectile 
center-of-mass.  In  other  words,  both  the  plane  of  the  projectile  motion  and  rotational 
kinematics  changes  as  a  result  if  the  target  interaction.  Projectile  rebound  merely 
requires  a  description  in  terms  of  the  projectile  center-of-mass,  i.e.,  the  plane  of  the 
projectile  motion  changes  as  a  result  of  the  target  interaction.  Bullet  ricochet  and  slug 
rebound  are  discussed  in  Sections  2.7.1  and  2.7.2,  respectively. 
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2.7.1  Bullet  Ricochet 


Bullet  ricochet  occurs  when  the  bullet  does  not  have  enough  incident  energy  to 
penetrate  the  target.  This  threshold  energy,  ,  is  dependent  on  the  projectile  and 

target  materials  as  well  as  their  shape  properties,  and  can  be  established  from 
experimental  results.  For  a  given  projectile-target  pair,  this  threshold  energy  can  be 
converted  to  a  threshold  impact  velocity,  ,  below  which  no  penetration  is  possible, 
and  the  bullet  ricochets  after  impact.  In  MRSAW  Vcr  has  been  calculated  from  Vcr  =  Vm 


such  that  the  exit  velocity  \s  Vex  =  0  in  Equation  (14).  This  gives:  = 


The  ensuing  motion  associated  with  bullet  ricochet  is  much  more  complicated  than  the 
initial  projectile  motion  prior  to  target  impact  when  no  bullet  ricochet  occurs.  Here  we 
face  another  aspect  of  the  bullet  motion  which  is  the  rotation  of  the  bullet  about  its 
center  of  mass,  G.  For  a  general  slender  body  (such  as  a  bullet)  the  kinematics  and 
kinetics  of  the  bullet  during  its  ricochet  is  presented  in  another  MRC  reference.® 

2.7.2  Slug  Rebound  Mode 

This  problem  is  analytically  identical  to  the  first  phase  of  the  projectile  exterior 
ballistics,  except  that  the  initial  conditions  are  adjusted  according  to  the  post  impact 
condition  of  the  bullet.  The  relevant  post  impact  conditions  are:  (1)  Rebound  speed 
after  impact;  and,  (2)  Direction  of  rebound  after  impact.  When  these  two  initial 
conditions  are  calculated  from  the  impact  equations,  equations  similar  to  Equations  (4) 
can  be  used  to  calculate  the  trajectory  of  the  projectile  after  impact. 

2.7.3  Coordinate  Systems  used  in  the  Input/Output  of  MRSAW 

All  coordinates  used  in  MRSAW  are  referred  to  the  initial  coordinate  system  used  in 
Phase  I  of  the  exterior  ballistics  portion  of  the  code  and  is  shown  in  Figure  1 1 .  All 
output  is  in  the  same  coordinate  system  for  first  and  subsequent  phases  of  projectile 
exterior  ballistics. 


2.7.4  Graphic  Representation  of  the  Results  from  MRSAW:  Input  and  Output  Files 

A  detailed  description  of  the  input  and  output  files  used  in  MRSAW  is  shown  in  Table 
13.  Descriptions  of  the  parameters  used  in  the  input  file  given  in  the  input  file  as 


^  Eisler,  R.D.,  and  Chatterjee,  A.K.,  Algorithm  Development  to  Describe  Retardation  in  Human  Tissue  of 
a  19.6  Grain  Flechette,  Mission  Research  Corporation  Report  MRC-COM-R-93-0372,  US  Army  Natick 
Technical  Report,  March  1993  (Unclassified  -  Limited  Distribution). 
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comments  to  the  right  of  the  numerical  input  along  with  the  units  used  in  the  code.  A 
typical  input  file  is  shown  below  in  Table  12. 


2.7.5  Buckshot:  Area  of  damage 

In  the  case  of  buckshot,  a  prescribed  number  of  spherical  projectiles  are  launched  with 
the  same  initial  muzzle  velocities  but  with  some  velocity  components  in  the  z-direction 
(Figure  1 1).  As  they  travel,  these  spherical  projectiles  cover  areas  proportional  to  the 
horizontal  distance  traveled.  The  constants  of  proportionality  for  particular  buckshot  is 
included  in  the  MRSAW  database,  and  can  be  used  to  determine  damaged  area  of  the 
target  for  this  type  of  ammunition.  For  the  second  phase  of  the  exterior  ballistics,  these 
results  are  shown  on  the  screen  when  relevant  during  the  execution  of  MRSAW.  The 
impact  locations  on  the  target  are  specified  in  MRSAW  in  cylindrical  coordinates.  The 
specific  radial  and  angular  coordinates,  in  the  target  plane,  of  each  striking  projectile  is 
determined  by  a  random  number  which  is  generated  consistent  with  the  pellet 
distribution  function,  n(R,  r),  defined  in  Equation  (15),  Section  2.7.9. 


Table  12.  Typical  Input  file  for  MRSAW  code 


3000.0,1.0,1200.0 

30.0,10.0 

10.0,5.0 

0.30 

1.0, 1.0, 1.0, 1.0, 1.0 
0.5, 0.5 
1.0, 1.0 


INITIAL  SPEED(FT/S) ,ANGLE(DEG) , RANGE (FT) 

CROSSWIND  SPEED (FT/S)  IN  THE  HORIZ.  &  VERT.  DIR. 
NORMAL  TGT.  ROT.  (DEG)  WRT.  Z-AND-Y-  OF  TGT.  AXES 
COEF.  OF  RESTITUTION  FOR  BULLET-TARGET  IMPACT 
BULL.  DEN.,  MASS(GH),  RADIUS,  BODY&TAIL  LENGTH (FT) 
FRONT  AND  REAR  TIP  DISTANCE  FROM  CM  G  (FT) 

BULLET  MI  ABOUT  G  (GM*FT**2)  AND  DIST.  G  FROM  TIP 


2.7.6  Fragmenting  Munitions 

Fragments  generated  during  bomb  explosions  have  less  deterministic  structures  than 
the  types  of  projectiles  discussed  before.  Some  of  the  problems  stem  from  the  fact  that 
the  size  and  shapes  of  these  fragments  can  at  best  be  described  by  probabilistic 
distribution  functions.  If  a  particular  fragment  is  identified  by  ‘the  so  called  UD  ratio’, 
which  is  the  ratio  of  the  longitudinal  dimension,  L,  to  the  lateral  dimension,  D,  then  the 
distribution  function  of  this  parameter  is  available  in  the  literature.  This  distribution 
function  describes  the  probability  of  UD  to  lie  within  a  given  domain.  In  order  to  use 
this  information  to  access  the  damage  produced  by  these  bomb  fragments,  we  need  to 
establish  either  the  deterministic  or  probabilistic  distribution  of  the  projectile 
trajectories.  It  is  obvious  from  the  deterministic  analysis  presented  in  the  preceding 
sections  for  the  trajectories  of  bullets  and  slugs,  that  we  need  to  know  the  initial 
conditions  regarding  fragment  ejection  velocities  and  ejection  angles.  No  theoretical 
development  is  available  for  the  partition  of  explosive  energy  for  individual  fragments  in 


# 


# 


# 
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Table  13.  Input  and  Output  Files  in  MRSAW 

X:  HORIZONTAL  DISTANCE(FT) 

Y:  VERTICAL  DISTANCE(FT); 

T:  TIME(SEC) 

XD:  HORIZONTAL  VELOCITY(FPS) 

YD:  VERTICAL  VELOCITY(FPS) 

FILE 

DESCRIPTION 

AIRBULL7.IN 

INPUT  FILE.  PARAMETERS  ARE  DECRIBED  IN 

TY1EX.PLT 

COMMENTS 

X  VS.  T  FROM  LINEAR  ANALYSIS 

TY2EX.PLT 

Y  VS.  T  FROM  LINEAR  ANALYSIS 

Y1Y1DEX.PLT 

X  VS.  XD  FROM  LINEAR  ANALYSIS 

Y1Y2DEX.PLT 

X  VS.  YD  FROM  LINEAR  ANALYSIS 

Y1Y2EX.PLT 

X  VS.  Y  FROM  LINEAR  ANALYSIS 

Y1Y2.PLT 

X  VS.  Y  FROM  NONLINEAR  ANALYSIS 

Y1Y1D.PLT 

X  VS.  XD  FROM  NONLINEAR  ANALYSIS 

Y1Y2D.PLT 

X  VS.  YD  FROM  NONLINEAR  ANALYSIS 

BULLRICO.IN 

INPUT  FILE  FOR  RICOCHETING  MODE  FOR 

BXY1.PLT 

BULLRICO.FOR 

OUTPUT  FROM  BULLRICO  CODE:  X  VS.  T  (FRONT  TIP) 

BXY2.PLT 

OUTPUT  FROM  BULLRICO  CODE:  Y  VS.  T  (FRONT  TIP) 

BXYG.PLT 

OUTPUT  FROM  BULLRICO  CODE:  X  VS.  T  (CENTER  OF 

AIRBULL7.DBG 

MASS) 

OUTPUT  FROM  BULLRICO  CODE:  DEBUG  FILE 

FTIP.PLT 

FINAL  OUTPUT:  X,Y,Z  LOCI  OF  FRONT  TIP 

RTIP.PLT 

FINAL  OUTPUT:  X,Y,Z  LOCI  OF  REAR  TIP 

CMG.PLT 

FINAL  OUTPUT:  X,Y,Z  LOCI  OF  CENTER  OF  MASS  G 

XYZALL.PLT 

FINAL  OUTPUT:  X,Y,Z  LOCI  OF  FRONT,  REAR  TIP  AND  G 

TXYZ12.PLT 

FINAL  OUTPUT:  T,X,Y,Z  LOCI  OF  FRONT  AND  REAR  TIP 

TYl.PLT 

FINAL  OUTPUT:  T,X,Y,Z  LOCI  OF  FRONT  TIP 

TY2.PLT 

FINAL  OUTPUT:  T,X,  Y,Z  LOCI  OF  REAR  TIP 

Y1Y1D.PLT 

FINAL  OUTPUT:  X,XD  OF  FRONT  TIP 

Y1Y2D.PLT 

FINAL  OUTPUT:  X,YD  OF  FRONT  TIP 

Y1Y2.PLT 

FINAL  OUTPUT:  X,Y  OF  FRONT  TIP 

an  exploding  munition  and  is  best  known  empirically  from  test  data.  Due  to  these 
problems,  we  concentrated  on  describing  the  probability  that  a  given  area  in  space  will 
be  impacted  by  fragments  of  sizes  within  some  range,  and  the  joint  probability  that 
these  fragments  will  have  some  critical  energy  or  striking  velocity  to  cause  damage  to 
physical  objects  (humans,  building  etc.)  occupying  this  area.  The  database  used  in  this 


32- 


Mission  Reseaich  Corporation 


M67004-96-C-0013  (2AB) 
MRC-COM-R-96-0461 


approach  is  existing  arena  test  data  published  by  the  DoD  in  Joint  Munitions 
Effectiveness  Manuals  (JMEM). 

2.7.7  Threshold  Velocities  for  Given  target/projectiie  pair 


For  specified  target-projectile  combinations,  it  remains  to  establish  critical  impact 
velocities,  ,  for  target  penetration.  We  also  need  to  validate  the  current  formulation 


for  v^,\.e.,  V, 


2Rd 


.  This  will  be  established  from  experimental  data  developed  in 


Phase  II.  The  Phase  II  database  will  be  incorporated  in  the  Phase  II  software  for 
accurate  determination  of  this  mode  of  motion  (rebound  or  penetration)  for  the  second 


phase  of  the  exterior  ballistics.  This  database  should  encompass  all  projectile-targets 
pairs  of  interest.  Currently,  MRSAW  requires  user  input  for  the  determination  of  the 
mode  of  the  second  phase  motion  both  for  bullets  and  slugs.  In  the  second  phase  of 


the  code  development  this  will  be  automated. 


2.7.8  Coefficient  of  Restitution 


In  the  case  of  a  projectile  rebounding  from  the  target  after  impact,  subsequent  motion 
of  the  projectile  depends  on  the  impulsive  force  and  its  direction  imparted  on  the 
projectile  during  impact.  The  time  integral  of  this  impulsive  force  over  the  time  of  impact 
is  the  impulse  which  is  equal  to  the  change  in  linear  momentum  (product  of  mass  and 
velocity).  This  is  a  vector  equation  that  relates  changes  in  components  of  momentum 
when  the  impulsive  force  over  time  of  impact  is  exactly  known.  However,  for  most 
problems  involving  impact,  this  information  is  not  easy  to  obtain  even  from  very 
controlled  experiments.  Since  frictional  forces  (or  lateral  force)  is  usually  small 
compared  to  the  force  exerted  to  the  projectile  in  the  direction  normal  to  the  target 
surface  at  the  point  of  impact,  the  momentum  equation  gives  a  conservation  of 
momentum  in  the  direction  tangential  to  the  surface  of  the  target.  To  calculate  the 
speed  and  direction  of  the  projectile  after  impact,  we  use  the  classical  law  of  impact 
which  states  that  the  velocity  of  separation  after  impact  is  equal  to  the  product  of  a 
constant  e  (called  the  coefficient  of  restitution)  and  the  velocity  of  approach  prior  to 
impact.  Experiments  on  impacting  bodies  show  that  e  depends  on  the  projectile-target 
pair.  In  the  development  of  MRSAW,  we  have  assumed  that  e  is  known  for  a  given 
projectile  target  pair.  In  Phase  II  we  will  experimentally  create  a  database  for  e 
employing  projectile-target  combinations  of  interest.  For  any  given  projectile-target  pair, 
e  can  be  calculated  by  measuring  the  speed  of  rebound  and  the  angle  of  rebound  of 
the  projectile  and  using  the  impact  law  and  momentum  equation  mentioned  above. 

2.7.9  Determination  of  Spray  Angle  and  Number  of  Pellets  within  a  given  radius  r 
at  a  given  range  R  from  Experimental  Data 

The  mathematical  formula  that  predicts  the  number  of  pellets  n  =n(R,r)  within  a  certain 
radius  rat  range  R  for  a  given  shot  in  a  burst  mode,  has  the  form: 
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n(R,r)  =  N.e 


(15) 


where  «and  ^are  model  parameters,  Rf  is  the  reference  range  and  N  is  the  total 
number  of  pellets  in  the  shot. 

This  formula  is  consistent  with  the  following  physical  conditions  that  need  to  be 
satisfied  by  the  function  n(R,r): 

1 .  At  very  close  range,  R  approaches  zero  and  n(R,r)  should  approach  N 

2.  At  a  given  range  R,  n(R,r)  should  be  an  increasing  function  of  the  radial  variable,  r. 

3.  r  should  lie  between  0  and  its  maximum  R  tan0  where  ^is  the  spray  angle. 

4.  For  any  given  range  R,  n(R,R  tanO)=N.  This  gives  fi^tanO 

Equation  15  for  n(R,r)  satisfies  all  these  conditions,  and  hence  is  a  candidate  function 
for  predicting  the  number  of  pellets  at  any  given  range,  R,  and  within  any  given  radius, 
r.  Its  applicability  depends  on  its  ability  to  predict  experimental  results  and  predicted 
trends  should  be  physically  plausible. 

To  apply  the  above  results  for  any  given  shot  with  known  number  of  pellets  N,  we  need 
to  have  at  least  two  experimental  results  where  the  number  of  pellets  n  within  a  given 
radius,  r,  at  a  given  range,  R,  are  known.  Let  these  two  data  points  be:  (1)  n=n1,  R=R1, 
r=r1;  and,  (2)  n=n2,  R=R2,  r=r2.  (The  experimental  results  below  were  obtained  from 
the  ammunition  manufacturer  and  were  reported  in  the  Sixth  progress  report  for  this 
effort). 

First  Experiment 

Shot  Size  1:  Pellet  Count  =  16:  Muzzle  Velocity  =1213  fps 
Energy  per  pellet  127.8  ft.-lbs 

Experimental  Result:  At  40  yards  range, 

4  pellets  within  a  radius  of  20  in. 

9  pellets  within  a  radius  of  30  in. 

Thus,  n1=4,  r1=20,  R1  =40x3x1 2  and  n2=9,  r2=30,  R2=R1 

N  :=  16 
n1  =  4 
r1  =  20 
R1  :  40- 3- 12 
n2  9 
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r2  :=  30 
R2  R1 
Rf  :=  40  3  12 
Rf  =  1.4410^ 


P  - 


r1 


'  n2 
Ini  — 

n; 


j2 

R2 


In 


In 


/n2\ 

ini 


P  =  0.026 
0  :=  atan(p) 

0  =  0.026 

.  I  N 

ln(  — 
n1 


a  := 


/R1 

Rf 


P  - 


r1 


/n1 

N 


a  =  116.774 


n(R.r) 


N  exp 


a- 


Rf 


Check  Results;  Theoretical  Prediction 
Check  Data: 


R  ^  40-3  12 
Rp  =37.095 
r  :=  20 
n(R,r)  =4 
r  :=  30 
n(R,r)  =9 

OK 
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Predictions: 

R  50-3-12 
R-p  =46.369 
r  :=  20 

n(R,r)  =  1.886 
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r  :=  30 

n(R,r)  =4.243 
r  :=  40 

n{R,r)  =9.546 


Second  Experiment: 

Shot  Size  4:  Pellet  Count  27;  Muzzle  Velocity  1268  fps 
Energy  per  pellet  67.3  ft.-lbs 

Experimental  Result;  At  40  yd.  range, 

9  pellets  within  a  radius  of  20  in. 

15  pellets  within  a  radius  of  30  in. 


N1  =  9 
N  =  27 
r1  =  20 
R  40-3-12 
N2  15 

R  =  1.44-10^ 
r2  =  30 
Rf  40-3-12 


Rf=  1.44-10 

r1  ,  /N2 

—  Ini  — 

R  \  N 


P  = 


t2 

R 


In 


N 


Inl 


/N2\ 


\N1/ 


p  =  0.029 
0  :=  atan(P) 

9  =  0.029 

Ji 

.N1 


Ini 


a 


R\l.  rl\ 

Rf  ^  R 


a  =  73.559 
n(R,r)  ;=  N-exp 


R  i  n 
a--— -(P  - 

Rf 
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Check  Results;  Theoretical  Predictions 


R  :=  40- 3- 12 
R-P  =41.507 
r  :=  20 
n(R,r)  =9 
r  :=  30 

n(R,r)  =  15  OK 


Other  Predictions:  At  a  range  of  50  yards, 


R  :=  50  3  12 
Rp  =  51.883 
r  :=  20 

n(R,r)  =5.297 
r  :=  30 

n(R,r)  =8.828 
r  :=  40 

n(R,r)  =  14.714 
r  45 

n{R,r)  =  18.996 
R  =  5-12 
R  =  60 
R-p  =  1.729 
r  :=  1.5 

n(R.r)  =26.685 

Plot  of  Pellet  Distribution  as  a  Function  of  r  for  any  given  range  R 

The  functional  form  of  the  distribution  of  pellet  within  a  circle  of  radius  r  at  range  R  is 
given  by: 


n{R,r)  =  N.e 


where  a  is  a  model  parameter  and  p=tane  where  0  is  the  spray  angle.  Rf  is  the  range 
where  experimental  data  is  collected.  In  this  case  Rf=40  yards=1440  inches. 


For  the  first  experiment  we  have  a=1 16.774,  p-0.026,  N-16 
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3  :=  0.026 
R  :=  40- 3- 12 

R  =  1.44-10^ 

N  :=  16 

Rf  =  1.44.10^ 

R.p  =37.44 
r  :=  0.75,1.50..  37.44 
R 


n1(r)  N  exp 


-a- 


Rf 


P  -  ^ 


For  the  second  experiment  we  have  a=73.559,  3=0.029,  N=27 


a  :=  73.559 
3  :=  0.029 
R  :=  40  3  12 

R  =  1.44-10^ 
N  :=  27 

Rf  =  1.44-10^ 


R-3  =41.76 
r  :=  0.72,1.44..  41.76 


n2{r) 


N-exp 


iP 


Graphical  Representation  of  Target  Damage  due  to; 

1 .  Burst  Mode  Shot  (3  rounds) 

2.  Automatic  Mode  Shots 

3.  Shots  with  Circular  Error  Probability  (CEP) 

2.7.10  Burst  Mode  Shots 

In  this  case,  the  distribution  of  number  pellets  n(R,r)  inside  a  circular  region  of  radius 
rat  any  given  range  R,  is  given  by  the  function: 


AllH 

n{R,r)  =  N.e  ^ 

For  any  given  shot,  the  model  parameters  a  and  /?  can  be  calculated  from 
experimental  data  on  the  number  of  pellets  for  at  least  two  pairs  of  values  of  R  and  r. 
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The  methodology  of  calculating  these  model  parameters  have  been  explained  in  the 
previous  sections.  Recall  that  p  is  related  to  the  spray  angle  0  through  p  =  \zxi.0  and 
N  is  the  total  number  of  pellets  in  the  shot. 

For  the  graphical  representation  of  damage  area  on  a  plane  target  located  at  a  range  R 
from  the  shot  location,  we  first  calculate  the  radius  of  coverage  from  =  i?tan^. 
Since  this  area  is  hit  by  N  pellets  at  random  locations,  we  first  generate  N  random 
numbers  within  0  and  1  using  a  random  number  generator.  The  i-th  number  (i=1,2,...,N) 

is  scaled  such  that  it  lies  between  and  r,  (ro=0)  where  r-  =  R  /?-ln 

is  the  inversion  of  Equation  (15)  and  Rf  is  the  range  used  in  the  experimental 
determination  of  a  and  /?from  Equation  (15).  These  numbers  give  the  radial 
coordinates  of  the  shot  points  inside  the  target.  To  determine  the  circumferential 
locations  for  these  shots,  we  generate  another  set  of  N  random  numbers  between  0 
and  2n:.  These  give  the  angular  coordinates  of  the  shots  on  the  target.  These  results 
are  then  shown  graphically. 


30 
27 
24 
21 
18 

n1(r 

-  15 

n2(r 

12 
9 
6 
3 
0 

Figure  12.  Number  of  Pellets  versus  Radius  at  a  Range  of  40  yards 
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2.7.11  Automatic  Mode 

In  the  case  of  automatic  mode,  we  assume  that  the  muzzle  velocities  remain  the  same 
and  only  the  shot  locations  may  vary  within  certain  range  while  the  shot  angles  varies 
between  some  maximum  and  minimum  angles  with  the  horizontal.  In  this  case,  the 
shots  on  target  should  lie  within  a  rectangular  region  where  the  vertical  boundaries  can 
be  established  from  the  maximum  and  minimum  shot  angles  from  locations  lying  on  the 
vertical  boundaries  on  the  shot  plane  using  the  analytical  code  MRSAW.  It  may  be 
recalled  here  that  for  the  determination  of  the  shot  location  on  a  given  target,  MRSAW 
requires  the  muzzle  velocity,  shot  location  and  shot  angle  with  the  horizontal. 

Horizontal  boundaries  will  be  identical  to  the  horizontal  boundaries  within  which  the 
automatic  shots  are  fired.  Once  the  boundaries  of  the  rectangular  shot  area  on  the 
target  are  determined,  graphical  representation  of  the  damage  area  on  the  target 
requires  the  generation  of  N  (N  is  the  number  of  shots  in  the  automatic  mode)  pairs  of 
random  number  generations  which  are  scaled  to  lie  within  the  rectangular  boundary. 
This  methodology  requires  only  two  runs  of  the  MRSAW  code  instead  of  the  N  number 
of  shots  actually  fired  from  the  shot  plane. 

2.7.12  Soldier  Stress  and  Aiming  Errors 

In  this  case,  the  damage  region  on  a  target  at  a  given  range  is  influenced  by  the  person 
firing  the  shots.  The  dimensions  describing  the  damage  depends  on  whether  the 
person  firing  these  shots  are  under  stress  or  not.  The  distribution  functions  describing 
the  linear  dimensions  of  the  damage  region  at  a  given  range  are  available  in  the 
literature  for  both  stressed  and  unstressed  people.  For  the  graphical  representation  of 
the  damage  zone,  we  need  to  generate  N  (number  of  shots  fired)  pairs  of  random 
numbers  lying  with  the  boundary  dictated  by  the  known  distribution  functions  at  a  given 
range.  MRSAW  has  been  programmed  to  read  this  data. 
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3.  VISUALIZATION  AND  SCENE  ANIMATION 

The  visualization  and  scene  animation  task,  Task  6.  identified  requirements,  tested 
methodologies,  and  evaluated  third-party  software  for  visualization  and  implementation 
of  effects  in  the  Phase  M  virtual  environment.  Two  strategies  were  considered  for 
developing  the  Phase  II  virtual  environment. 

The  first  approach  was  to  develop  wrap-around  code  for  an  existing  commercial 
visualization  system.  Several  candidate  commercial  systems  were  considered  for  this 
application  including  MultiGen  ™,  Coryphaeus  ™,  and  Division  ™. 

The  second  approach  was  to  custom  design,  develop,  and  integrate  a  three 
dimensional  modeling  and  real-time  simulation  system.  This  system  would  be  tailored 
specifically  for  the  types  of  objects,  environments,  and  scenarios  which  would  be 
simulated  in  the  Phase  II  system.  Both  approaches  would  yield  compelling  simulations 
of  the  virtual  environment  however  implementation  of  this  second  approach  was  quickly 
seen  to  exceed  resources  available  in  the  Phase  II  effort. 

The  proposed  approach  is  a  hybrid  of  the  above  strategies  and  develops  the  virtual 
environment  and  interfaces  by  integrating  commercially  available  off-the-shelf  hardware 
and  software.  During  Phase  I,  representative  hardware  and  software  was  evaluated  for 
compatibility  with  the  effects  to  be  implemented  in  the  Phase  II  virtual  environment.  The 
effects  to  be  implemented  were  flowed  down  from  the  Scenario  Development  task.  Task 
1,  the  Damage-Response  Mode  Matrix  task.  Task  2,  the  Weapons  Effects  task.  Task  3, 
and  the  Structural  Integrity  task.  Task  4. 

3.1  SOFTWARE  EVALUATION 

Existing  commercial  software  systems  were  evaluated  to  understand  constraints 
associated  with  the  real-time  system  development  in  Phase  II  and  to  gain  insight  as  to 
how  current  state-of-the-art  technology  in  real  time  modeling  and  simulation  could  be 
exploited  for  the  Phase  II  effort.  This  would  enable  us  to  determine  whether  off-the- 
shelf  software  could  be  used  for  this  project  and  also  what  deficiencies,  if  any,  exist  in 
the  software.  While  many  real  time  modeling  and  simulation  packages  are  available, 
we  focused  our  evaluation  on  two  competing  packages,  MultiGen  Inc.  (San  Jose,  CA) 
and  Coryphaeus  Software  Inc.  (Los  Gatos,  CA).  Both  systems  were  chosen  because  of 
their  widespread  use  within  the  Department  of  Defense. 

The  main  criteria  used  to  determine  the  effectiveness  of  the  software  (in  no  particular 
order)  were; 

•  modeling  features 

•  texture  mapping 


M67004-96-C-0013  (2AB) 
MRC-COM-R-96-0461 


Mission  Research  Corporation 


•  dynamic  lighting  capabilities 

•  special  effects  (explosions,  smoke,  fire,  debris,  flak,  etc.) 

•  animation  features 

•  audio  capabilities 

•  DIS  compatibility 

•  ease  of  learning  and  ease  of  use 

•  programmability  and  customization 

•  support  for  other  data  formats 

•  real  time  data  from  external  application 

•  portability 

3.1.1  MultiGen™ 

The  first  tool  evaluated  was  MultiGen.  While  Multigen  proved  to  be  a  very  strong  3D 
modeling  package  for  real-time  simulations  and  satisfied  many  of  the  established 
criteria.  The  special  effects  available  as  options  in  the  software  were  limited  for  our 
application.  Once  3D  models  were  completed,  there  was  no  easy  way  to  view  them  in 
real-time  due  to  a  lack  of  a  simulation  development  environment.  MultiGen  is  primarily 
a  real  time  3D  modeling  tool  and  not  a  real  time  simulation  development  tool.  It  was 
determined  that  we  needed  to  evaluate  real  time  simulation  development  environments 
in  addition  to  3D  modeling  environments. 

3.1.2  Coryphaeus  ™ 

Designer’s  Workbench™  (3D  modeling  package)  and  EasyScene™  (real  time 
development  environment)  from  Coryphaeus  Software  Inc.  were  evaluated.  Designer’s 
Workbench  proved  to  be  a  good  modeling  package  for  real-time  databases.  While 
somewhat  limited  in  complex  geometry  creation  tools,  it  has  sufficient  tools  for  creating 
the  simple  geometry  required  for  real-time  applications.  Also,  the  ability  to  read 
geometry  from  a  variety  of  other  modeling  packages  and  other  real-time  applications 
helps  make  up  for  modeling  deficiencies. 

Texture  mapping  is  straight  forward  using  the  graphical  user  interface.  Unfortunately, 
even  a  texture  mapped  cube  requires  more  graphics  power  than  is  available  on  a  SGI 
lndigo2  Extreme™  which  we  were  using  during  the  Phase  1  activity.  This  proved  to  be 
the  limiting  factor  in  evaluating  texture  mapping  performance. 

Most  of  the  lighting  options  are  geared  more  toward  aircraft  simulations  but  it  is 
possible  to  simulate  lighting  conditions  inside  buildings  as  well.  To  add  depth  to  the 
simulation,  shadows  can  be  activated.  Infrared  lighting  conditions  can  also  be 
simulated  without  much  difficulty. 

As  with  the  other  real  time  3D  packages  tested.  Coryphaeus  has  a  limited  selection  of 
special  effects.  However,  additional  special  effects  modules  can  be  written  in  C/C++ 
and  incorporated  into  the  simulation.  Since  our  application  requires  debris  and  ejecta 
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for  a  wide  variety  of  materials  and  munitions,  “canned”  special  effects  would  not  be 
utilized  exclusively  anyway. 

Audio  effects  can  be  attached  to  objects  or  events  and  triggered  within  a  simulation. 
Thus,  it  will  be  possible  to  record  the  sounds  from  actual  explosions  and  play  them 
during  the  triggered  event.  This  will  enhance  the  perceived  realism  of  the  experience. 

Coryphaeus  claims  to  have  built  in  DIS  capability  to  allow  access  to  a  DIS  network 
without  writing  additional  code.  DIS/HLA  requirements  are  undergoing  massive 
revisions  however  and  these  issues  will  have  to  be  re-examined  during  Phase  II  as  new 
requirements  unfold  and  as  the  wrap-around  software  is  being  developed. 


Programming  additional  features  and  creating  user  defined  modules  will  be  a  very  key 
part  of  the  phase  II  effort.  Coryphaeus  has  hundreds  of  user  functions  which  can  be 
called  from  a  C/C++  program  and  incorporated  into  the  simulation.  Coryphaeus  is  also 
totally  compatible  with  IRIS  Performer  ™  (Silicon  Graphics  ANSI  C  and  C++ 
programming  interface  for  creating  real-time  visual  simulation  and  other  interactive 
graphics  applications)  which  offers  additional  functions  and  features.  There  is  also 
support  for  multiple  channels  of  input/output  data,  multiple  CPU’s,  and  multiple 
graphics  pipes  (to  feed  graphics  to  multiple  display  devices).  These  features  will  enable 
optimization  of  the  hardware  and  software  architecture  environment  for  real  time 
performance. 

One  very  strong  feature  of  Coryphaeus  is  the  ability  to  import  data  from  other  formats 
and  applications.  This  enables  us  to  take  advantage  of  the  wide  variety  of  3D  geometry 
which  have  been  created  for  other  applications.  Geometry  can  be  imported  from  about 
30  other  modeling  and  real-time  applications  including  Wavefront  ™  (advanced 
visualization  and  inverse  kinematics  software  commonly  used  in  scientific 
visualization).  Alias  ™  (advanced  surfacing  and  animation  software  commonly  used  in 
design  and  entertainment),  IRIS  Performer  ™  (Silicon  Graphics  ANSI  C  and  C++ 
programming  interface  for  creating  real-time  visual  simulation  and  other  interactive 
graphics  applications),  Lightscape  ™.  (A  visualization  system  combining  proprietary 
radiosity  algorithms  with  a  physically  based  lighting  interface  commonly  used  in  the 
architecture  industry).  The  extensive  array  of  3D  models  available  from  Image  data  can 
be  imported  for  use  as  textures  and  special  effects.  Motion  and  event  data  can  also  be 
imported  and  applied  to  objects  in  the  database.  This  will  be  crucial  as  the  algorithms 
which  describe  weapons  effects  will  be  incorporated  in  this  way.  The  only  computer 
platform  Coryphaeus  software  currently  runs  on  is  Silicon  Graphics  but  this  is  a 
limitation  shared  by  many  other  visualization  tools  as  well. 

3.2  DEVELOPMENT  APPROACH 

Several  o/f-f/7e-s/?e/f  solutions  exist  for  creation  of  3D  geometry  for  real  time 
simulations.  Both  Coryphaeus  and  Multigen  were  capable  of  creating  the  objects 
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necessary  for  MOUT  type  virtual  environments.  It  was 
discovered  however  that  the  modeling  activity  is  a  less 
significant  element  of  the  overall  simulation  than  initially 
perceived.  That  is,  development  of  the  weapons  effects 
simulation  is  much  more  complex  and  represent  more 
stressing  constraints  on  Phase  II  development. 


Figure  13.  Ceiling  The  software  evaluated  contained  tools  for  creating 

Surface  following  Ballistic  (aircraft  simulations),  weather  conditions 

Test  with  12-guage  vjgjbjijty  background  color),  and 

shotgun  using  Number  4  triggerable  effects  (explosions,  flak,  and  debris), 

buckshot  from  8  feet.  Explosions  are  represented  as  reddish  orange  plumes  of 

flames  and  smoke  that  dissipate  from  a  model.  Flak  is 
represented  as  a  light  gray  and  black  plume  that  dissipates  from  the  model.  Debris  is 
viewed  as  small  particles  that  move  outward  from  the  origin  of  the  model.  While  these 
triggerable  effects  are  adequate  when  viewed  from  the  far-field  (from  an  airplane  to  the 
ground  for  example),  they  are  not  as  useful  for  near  field  views  of  weapon  effects 
prevalent  in  an  urban  warfare  application  (refer  to  Figure  13  and  Figure  14). 


To  enhance  the  realism  of  the  effects,  additional  programming  will  be  necessary.  We 
can  build  on  the  base  effects  and  add  the  functionality  necessary  to  create  convincing 
simulations  of  our  weapons  effects.  This  will  be  accomplished  by  developing  ANSI  C 
and  C++  routines  which  call  functions  from  Coryphaeus  API,  IRIS  Performer,  X,  and 
OpenGL  ™  (platform  independent  graphics  libraries  and  functions  developed  by  SGI) 
executable  libraries.  These  routines  are  then  compiled  and  attached  to  the  simulation 
at  run-time. 


In  conclusion,  a  hybrid  approach  that  uses  off-the-shelf  commercial  software  in 
combination  with  custom  developed  software  will  be  necessary.  The  Coryphaeus 
products  evaluated  will  form  the  foundation  for  the  development  environment.  Most,  if 
not  all,  of  the  3D  environment  will  be  created  using  Designer’s  Workbench.  EasyScene 
will  be  used  to  create  the  basic  simulation  environment.  Additional  programming  will  be 


Figure  14.  Backface  Debris  and  Ejecta  from  Interior  Sheet  Rock  Wall  subject  to  12 
gauge  shotgun  with  #4  Buckshot  fired  from  8  feet 


(STRICOM  Phase  I  Ballisitc  Experiments  conducted  by  MRC) 
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implemented  to  develop  the  weapons  effects  required  by  urban  warfare  and  the 
methodology  necessary  to  link  the  weapons  effects  algorithms  with  the  weapons  effects 
visuals. 

3.3  McKenna  MOUT  SITE  DEVELOPMENT 

For  purpose  of  software  evaluation,  the  McKenna  MOUT  Town  Hall  building  (shown  in 
Figure  15)  was  created  using  Designer’s  Workbench  from  Coryphaeus.  The  geometry 
accurately  represents  the  engineering  drawings  as  supplied  by  the  United  States  Army 
Infantry  School  at  Fort  Banning.  All  walls,  windows,  and  doors  on  all  three  floors  have 
been  accurately  positioned  and  the  stairwell  in  the  center  of  the  building  allows  access 
to  all  floors.  The  floor  plan  for  the  ground  floor  is  illustrated  in  Figure  16.  Figure  17 


Figure  15.  Exterior  of  McKenna  MOUT  Site  Town  Hall  Building 


illustrates  the  main  entrance  of  the  McKenna  MOUT  Town  Hall  Building,  and  Figure  18 
shows  the  stairwell. 


With  some  effort,  the  entire  building  can  be  explored  in  near  real-time.  A  simple  fly- 
through  of  the  McKenna  MOUT  Town  Hall  Building  was  provided  to  STRICOM  in  VHS 
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format.  The  original  intent,  as  is  the  case  with  all  real-time  efforts,  was  to  utilize  texture 
mapping  as  much  as  possible  to  minimize  the  need  for  geometry.  However,  the 
graphics  requirements  of  texture  mapping  exceeded  the  performance  capabilities  of  the 
SGI  Indigo^  Extreme  used  during  the  Phase  I  effort.  Even  a  cube  with  a  brick  texture 
map  rendered  the  machine  nearly  useless  for  development.  The  cube  rotated  very 


Figure  16.  First  Floor  McKenna  Town  Hall  Layout 


slowly  and  was  jerky  and  all  menus  responded  several  seconds  after  selection.  Thus, 
more  geometry  was  created  and  simple  materials  were  used  to  represent  surfaces.  No 
textures  were  used.  This  is  only  a  temporary  problem  as  any  implementation  of  the  final 
system  would  undoubtedly  run  on  a  high  performance  computer  system  (SGI  High 
impact  or  better).  When  interior  photos  are  available  and  satisfactory  hardware  is 
available,  a  very  realistic  representation  of  the  McKenna  MOUT  site  is  achievable.  The 
Phase  I  building  is  meant  to  be  used  as  a  representative  MOUT  object  not  an  exact 
replica. 
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To  conserve  resources  during  Phase  I,  several  issues  necessary  for  a  compelling 
weapons  effects  simulation  were  investigated  decoupled  from  the  algorithms  describing 
physical  effects  and  the  Town  Hall  fly-through  developed  in  this  and  other  tasks  of  the 
Phase  I  effort.  That  is,  the  weapons  effects  simulations  were  evaluated  for  availability 
and  effectiveness  but  not  developed  to  the  point  where  the  simulations  could  be 
integrated  into  the  Town  Hall  fly-through.  This  was  done  to  highlight  requirements  that 
need  to  be  imposed  on  the  Phase  II  virtual  environment  and  to  allow  adequate  time  to 


Figure  17.  Town  Hall  Main  Entrance 


evaluate  off-the-shelf  software  in  meeting  these  requirements.  Dynamic  lighting, 
animation  effects  (opening  and  closing  doors,  etc.),  collision  detection,  and  level  of 
detail  were  evaluated  in  this  manner.  For  example,  a  simple  model  of  a  door  and  a  wall 
was  created  to  investigate  the  mechanism  of  animating  the  opening  and  closing  of  the 
door.  Effects  or  functions  requiring  texture  mapping  (smoke,  dust,  fire,  fabric,  wall 
hangings,  etc.)  were  not  evaluated  due  to  the  hardware  limitation  previously  mentioned. 

3.4  Weapons  Effects  Visualization 

The  primary  elements  of  convincing  weapons  effects  in  an  immersed  environment  are; 
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•  sound  effects  (of  weapon  fire  and  impact  with  objects) 

•  dust,  flak  (airborne  particulate  created  at  bullet  entry  and  exit  from  object) 

•  debris  (pieces  of  objects  remaining  after  impact  either  around  penetration  hole  or  on 
the  floor  below  the  hole) 

•  smoke  and  fire  (when  appropriate) 

•  visible  change  of  geometry  (actual  holes  in  objects  or  surfaces 


Figure  18.  Town  Hail  Stairwell 


If  any  of  these  factors  are  missing,  the  simulation  will  be  less  effective.  Sound  effects  of 
weapons  triggered  at  the  appropriate  time  are  not  a  problem.  Audio  recordings  from 
ballistic  experiments  will  be  adequate  for  first  level  realism.  However,  all  other  factors 
indicated  above  will  require  additional  development  outside  of  the  Coryphaeus 
environment.  While  some  of  the  necessary  controls  are  available  for  dust,  smoke,  and 
fire,  they  were  not  considered  sufficiently  robust  to  handle  the  variety  of  munitions  and 
material  combinations  possible  in  the  MOUT  environment. 
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The  debris  resulting  from  an  explosion  or  from  a  weapon  being  fired  into  a  surface,  as 

•  well  as  the  visible  change  in  the  affected  surface  geometry  present  the  biggest 
challenge.  What  level  of  realism  is  necessary?  First  level  realism  can  be  achieved  by 
digitally  pasting  a  photo  of  an  actual  bullet  hole  onto  the  surface  where  the  bullet 
penetrated  (Figure  19). 

•  The  photo  (texture)  would  come  from  a  library  of  digitized  photos  obtained  from  ballistic 
experiments  of  various  materials  and  munitions.  The  weapons  effect  would  be 
completed  by  adding  the  effect  of  dust,  debris,  or  ejecta  resulting  from  the  bullet. 


Figure  19.  Texture  Mapped  Surface  of  Impact  Surface  of  Interior  Wall  with  12 
gauge  Shot  Gun  Hole  (Number  4  buckshot  and  range  of  8  feet). 


A  problem  with  this  approach  is  that  it  can’t  account  for  the  situation  where  sizable 
holes  are  created  by  the  weapon.  For  example,  if  a  hole  large  enough  to  look  through 
into  the  next  room  is  strategically  placed  in  a  wall,  a  texture  map  cannot  alone  be  used 
for  the  simulation.  Applying  the  texture  of  a  hole  to  a  surface  does  not  create  a  hole  in 
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the  surface.  In  that  way,  a  texture  map  can  be  thought  of  as  wallpaper.  If  the  wallpaper 
does  not  have  a  physical  hole,  placing  it  on  a  surface  will  not  create  a  hole. 

The  only  way  to  simulate  this  situation  and  achieve  a  higher  level  of  realism  is  by 
geometric  modification.  If  we  knew  where  to  hole  was  to  be  located  each  time  the 
simulation  was  entered,  we  could  create  a  hole  in  the  surface  in  advance  and  “uncover” 
it  at  the  appropriate  time.  Since  a  weapon  can  be  fired  at  any  surface  in  the  simulated 
environment  from  any  direction,  the  possible  locations  for  the  “hole”  are  innumerable. 
Thus,  we  must  develop  a  method  to  reconfigure  the  geometry  of  the  affected  surface  in 
real  time. 


Figure  20.  Interior  Wall  with  Hole  from  12  gauge  shotgun  and  Number  4  Buckshot 

Adding  the  texture  (photo)  of  the  hole  will  provide  the  additional  detail  necessary  for  the 
debris  around  the  hole  (refer  to  Figure  21).  The  addition  of  an  effect  to  simulate  dust, 
debris,  or  ejecta  from  the  penetration  will  complete  this  effect. 
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The  weapons  effects  algorithms  developed  in  Tasks  3  and  4  can  already  predict  the 
hole  pattern  resulting  from  various  weapons  and  material  combinations.  Also,  data  is 
available  from  the  simulation  to  precisely  locate  in  3-dimensional  space  where  the 
weapon  was  fired  from  and  the  coordinates  of  each  object  in  the  environment.  With  this 
information  available,  it  is  possible  to  develop  a  function  which  will  take  the  geometry  of 
a  wall  for  example,  and  re-create  the  wall  to  contain  the  hole  created  by  the  weapon. 
This  new  wall  geometry  would  then  be  passed  back  to  the  simulation  where  it  will 
replace  the  original  wall  until  the  simulation  is  re-initialized. 


The  wall  with  the  hole  (see  to  Figure  20)  has  been  created  by  digitally  “tracing”  the  hole 
in  the  photo  and  removing  that  portion  of  the  wall.  It  remains  to  be  determined  what 
level  of  the  geometric  reconfiguration  calculations  can  be  performed  in  real  time.  As  the 


Figure  21.  Object  with  Geometric  Hole  and  Applied  Texture  Map 


simulation  progresses,  more  and  more  weapons  will  be  fired,  resulting  in  more 
weapons  effects  being  simulated.  This  creates  a  potential  problem  for  real  time 
geometric  reconfiguration  calculations  and  real  time  display  of  the  resulting  3D 
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environment.  Figure  22,  Figure  23,  and  Figure  24  illustrates  the  accumulative 
polygonal  effect  created  by  adding  more  and  more  bullet  holes  to  an  object.  Figure  22 
is  the  polygonal  representation  of  wall.  Figure  23  is  the  polygonal  representation  of  a 
wall  with  one  bullet  hole,  and  Figure  24  is  the  final  wall  with  multiple  holes. 

Note  that  the  polygon  count  has  increased  from  6  to  807  polygons.  If  50  bullets  holes 
were  created  instead  of  5,  the  polygon  count  would  be  over  7000!  This  increase 
dramatically  illustrates  why  texture  mapping  is  an  attractive  option  for  far  field  viewing. 
Graphics  hardware  can  display  50  texture  mapped  bullet  holes  much  easier  than  an 
additional  7000  polygons.  It  is  very  important  to  use  geometric  reconfiguration 
cautiously.  The  algorithm  would  need  to  be  robust  enough  to  determine,  based  on 
viewing  distance,  whether  a  texture  map  should  be  used  (Figure  19)  instead  of  actual 
geometry  (Figure  21). 


i 

i 


Figure  22.  Solid  Wall  -  6  Polygons 


Another  option  is  to  use  simplified  bullet  holes.  For  the  bullet  hole  in  Figure  20  we 
“traced”  the  outline  of  the  bullet  hole.  We  could  instead  use  a  6-10  sided  circle  for  the 
hole  and  still  create  a  hole  in  the  surface.  This  is  not  as  visually  compelling  but  would 
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reduce  the  polygon  count  significantly.  This  option  could  be  used  as  the  viewer  walks 
toward  the  wall.  The  actual  hole  could  be  substituted  as  the  viewer  reaches  the  wall. 

If  the  initial  geometry  of  all  walls  and  objects  in  the  environment  is  kept  simple 
(minimum  curvature,  minimum  surface  detail)  and  a  dedicated  processor  is  used  for  the 
calculations  (this  can  even  be  a  PC  based  computer),  satisfactory  geometric 
reconfiguration  performance  can  be  achieved.  Also,  implementing  level  of  detail  control 
into  the  geometry  reconfiguration  algorithm  (realism  of  bullet  hole  shape  as  a  function 
of  viewer  distance  from  the  surface)  will  help  minimize  polygon  count  and  thus  increase 
display  performance. 
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Figure  23.  Wall  with  One  Bullet  Hole  -  159  Polygons 


3.5  Hardware  and  Software  Architecture 

The  overall  simulation  structure  can  be  designed  to  execute  as  a  single  UNIX  process. 
However  to  permit  real  time  performance,  an  environment  with  more  than  one 
processor  is  desirable.  A  high  end  SGI  computer  (SGI  High  Impact  or  better)  will  be 
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required  for  the  graphics  display  of  the  simulation.  To  achieve  real  time  performance,  it 
is  important  to  perform  as  many  auxiliary  computations  (weapons  effects  algorithms, 
geometric  reconfiguration  function,  etc.)  as  possible  on  separate  processors.  The 
additional  processors  could  be  Pentium  ™  based  PC’s.  As  more  characters  are 
immersed  in  the  simulation,  additional  graphics  processors  would  be  required. 

Each  immersed  character,  in  addition  to  having  some  type  of  tracking  device,  will  also 
have  one  or  more  weapons.  At  this  time,  no  significant  effort  has  been  expended  to 
evaluate  these  devices.  The  core  system  proposed  in  phase  one  has  built  in 
capabilities  to  handle  the  input  from  these  devices  and  any  specific  issues  will  be 
revisited  during  Phase  II. 


Figure  24.  Wall  with  Five  Bullet  Holes  -  807  Polygons 


In  the  attached  figure,  the  grayed  areas  indicate  the  focus  of  our  effort  and  extensions 
to  the  state-of-the  art.  The  non-grayed  areas  indicate  existing  off-the-shelf  solutions 
that  will  be  exploited  and  adapted  for  our  purposes.  Three  dimensional  visual 
representations  of  buildings  in  general  and  the  McKenna  MOUT  site  in  particular  are 
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plentiful.  However,  models  suitable  for  real-time  fly-through  visual  simulations  are  not 
as  readily  available.  A  baseline  model  of  the  McKenna  Town  Hall  building  was 
developed  in  Phase  I.  A  simulation  using  this  model  can  accommodate  simple 
animation  involving  simple  translation  and  rotation  of  objects.  The  entire  McKenna 
MOUT  site  can  be  developed  in  a  similar  manner  if  necessary. 

The  problem  for  the  effects  we  wish  to  simulate  are  threefold  in  terms  of  extensions  to 
the  state  of  the  art.  First,  we  need  to  allow  “atomic”  elements  of  the  simulation  to 
deform  in  a  non-rigid  manner  (e.g.,  fragment,  crush,  bend,  etc.).  Existing  visual 
simulations  can  not  do  this  since  they  merely  render  visible  surfaces.  Second  we  need 
to  employ  data  bases  with  fundamentally  different  spatial  and  temporal  resolutions  to 
accomplish  the  above.  Further,  many  of  these  effects  occur  in  parallel  and  must  be 
appropriately  queued  to  the  number  and  speed  of  the  hardware/processors  being  used. 

Finally,  various  decisions  must  be  made.  First,  what  do  we  need  to  show  in  order  to 
make  the  simulation  visually  compelling.  In  addition  to  visuals,  what  level  of  audio  or 
other  non-visual  cues  are  necessary.  Second,  what  don’t  we  have  to  show  but  never- 
the-less  need  to  know.  The  “atomic”  or  fundamental  elements  of  the  simulation  must  be 
also  be  identified.  This  last  issue  involves  various  trade-offs  in  terms  of  simulating  rigid- 
body  and  non-rigid  body  processes,  the  number  of  polygons  in  the  visual  simulation, 
and  real-time  performance.  Also  what  is  the  appropriate  level  of  approximation  and 
packaging  of  the  physical  processes  being  simulated. 

After  wrestling  with  these  problems,  it  became  apparent  that  the  weapons  effects 
models  cannot  be  separated  from  the  visualization  strategy.  The  visualization  of  these 
effects  is  the  real-time  constraint.  The  physical  effects  models  would  therefore  have  to 
be  tailored  both  in  terms  of  level  of  approximation  and  in  terms  of  model  parameters  so 
that  that  a  minimum  number  of  atomic  elements  and  polygons  occur  in  the  visualization 
so  that  the  visualization  of  weapon  effects  can  be  done  in  real-time. 

For  example,  in  determining  the  effect  of  a  shot  gun  blast  through  an  interior  wall.  The 
simulation  will  not  perform  the  detailed  calculations  that  describe  the  interaction  of  the 
projectiles  with  each  of  the  wail  materials.  Rather,  a  library  of  interior  wall  constructions 
will  be  provided  with  parametric  solutions  specific  to  the  projectile  and  wall  construction 
combination.  The  parametric  solutions  will  be  compared  with  detailed  calculations  and 
experiments  in  Phase  II.  For  some  scenarios  we  have  already  done  this  in  Phase  I. 

Similarly,  debris  fields  will  require  physical  models  that  describe  particulate  density, 
and  visual  transmittance  as  a  function  of  time.  For  debris  fields  with  casualty  producing 
debris  (e.g.,  glass  shards),  it  is  not  practical  (due  to  limitations  associated  with  the  3D 
visualization  methodology  of  debris  fields)  to  track  each  shard  but  rather  develop  a 
stochastic  model  which  describes  the  distribution  of  ejecta  and  ejecta  velocities,  and 
performs  a  random  draw  to  determine  striking  conditions  on  personnel.  Once  the 
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striking  conditions  are  prescribed,  wound  calculations  (to  the  level  of  fidelity  required 
for  non-medical  purposes)  can  be  done  in  real-time. 

We  already  have  a  mature  methodology  to  perform  wound  analysis  on  personnel  that 
operates  faster  than  real-time  and  we  merely  have  to  adapt  this.  It  is  likely  that  most  of 
the  software  we  will  need  will  be  developed  by  us  either  in  a  new  contract  from 
Natick/Simulation  Technologies  or  Natick’s  Phase  II  SBIR.  Since  we  already  know  what 
to  do  here  relative  to  penetrating  wounds  we  did  not  waste  Phase  I  effort  belaboring 
this  issue  but  we  will  describe  our  approach  in  the  Phase  II  proposal.  As  part  of  the 
STRICOM  effort  we  only  expect  to  focus  on  penetrating  wounds.  In  the  Natick  Phase  II 
effort  we  will  be  focusing  on  other  types  of  trauma  which  can  be  exported  to  the 
STRICOM  simulation.  For  medical  simulation  we  can  provide  interfaces  with  the 
DARPA  surgical  simulators  that  are  already  under  development  and  which  eventually 
will  be  operating  in  real-time  (there  is  no  point  in  reinventing  the  wheel).  In  other  words 
when  someone  approaches  a  casualty  and  looks  down  within  a  prescribed  range  (you 
would  switch  to  the  DARPA  simulators). 

3.6  STRICOM  Phase  II  SBIR  Software  Architecture 

Input  files  will  be  monolithic  (i.e.  containing  both  visual  and  structural  information)  to 
maintain  synchronization  of  the  data  as  the  number  and  size  of  the  data  bases  grow. 
We  expect  that  the  input  format  would  permit  interface  with  other  visualization  tools. 

The  overall  structure  will  be  designed  to  execute  as  a  single  UNIX  process.  However, 
in  an  environment  with  more  than  one  processor  available,  various  secondary 
processes  can  be  forked  to  permit  better  real  time  performance.  We  will  build  a 
hardware  environment  with  multiple  SGI  Indigos  and  a  SGI  High  End  Impact  for 
development  purposes.  Toward  the  end  of  the  effort  we  will  export  the  software  to  a 
STRICOM  or  Natick  hardware  environment  with  higher  end  platforms  to  verify  real-time 
performance  and  visual  fidelity  (this  will  require  someone  from  our  organization  to  be 
resident  at  Natick  and/or  STRICOM  for  up  to  two  months  at  the  end  of  the  effort. 

We  intend  to  provide  an  environment  with  multiple  characters  interacting  (initially  two, 
this  is  merely  limited  by  the  number  of  processors  and  input  devices  available).  Two 
characters  will  allow  us  to  confront  most  of  the  fundamental  problems  without  spending 
a  lot  of  funding  on  hardware. 

For  any  simulation  system  there  has  to  be  a  core  structure  which  maintains  the  real 
time  state  of  the  simulated  environment.  Using  EasyScene™  from  Coryphaeus™ 
Software  will  provide  a  foundation  on  which  to  build  the  “Steady  State  Simulation 
Process”.  EasyScene  will  provide  real  time  video  output  from  an  SGI  platform.  The 
Steady  State  Simulation  Process  controls  all  of  the  routine  simulation  effects  which  only 
require  translation  and  rotation  of  objects  (Rigid  Body  Motion)  contained  in  the  data 
base. 
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An  “event”  is  defined  to  be  any  occurrence  in  the  simulation  which  may  physically  alter 
one  or  more  “atomic  structures”  in  the  data  base  (Non-Rigid  Body  Deformation).  An 
atomic  structure  is  defined  to  be  one  object  or  group  in  the  data  base  which  has  a 
defined  structural  composition  (e.g.  a  window  or  a  wall  panel). 

When  the  Steady  State  Process  detects  that  an  event  has  occurred,  the  event  is 
queued  for  handling  by  an  “Event  Process.”  There  may  be  more  than  one  Event 
Process  available  to  handle  events,  depending  on  the  system  resources  available. 
Some  Event  Processes  will  require  a  dedicated  CPU  to  maintain  real  time  performance. 

The  function  of  the  Event  Process  is  to  evaluate  the  nature  of  the  event  and  the  data 
base  environment  in  which  it  occurred.  From  the  center  of  the  event,  the  environment  is 
evaluated  and  any  atomic  objects  which  may  be  affected  by  the  event  are  queued.  The 
interaction  of  the  object  and  the  event  not  only  affects  the  object,  but  may  alter 
subsequent  objects  that  are  being  checked.  This  permits  shadowing  effects  for  blasts 
and  the  determination  of  penetration  effects  of  firearms. 

The  Event  Process  determines  those  objects  which  will  be  affected  first  (the  objects 
closest  to  a  blast  or  adjacent  objects  in  the  line  of  fire)  and  queue  the  appropriate 
combinations  of  “Structural  Processes”  and  “Character  Processes.”  After  these 
processes  have  evaluated  the  effect  on  themselves  and  the  affect  of  the  event,  the  data 
is  returned  to  the  Event  Process.  If  the  event  can  still  affect  another  object,  the  next 
object  is  queued  for  evaluation.  This  continues  until  the  effect  of  the  event  dissipates. 

The  Structural  Process  has  two  functions,  determine  the  effect  of  an  event,  and  create 
a  new  visual  representation  of  the  object  as  it  would  appear  after  the  event.  The  visual 
representation  is  made  up  of  two  parts,  the  permanent  representation  and  any  transient 
effects  (e.g.  dust,  smoke,  debris).  The  updated  visual  representation  is  passed  back  to 
the  Steady  State  Process  for  rendering. 

The  base  requirements  for  the  Character  Process  are  the  same  as  the  Structural 
Process  with  the  addition  of  character  animation  and  biomechanical  functions.  These 
biomechanical  functions  would  be  emphasized  in  the  Natick  Phase  II  SBIR  and 
developed  to  a  much  higher  level  of  sophistication  and  fidelity  than  in  the  STRICOM 
effort.  Similarly,  the  STRICOM  effort  would  emphasize  the  structural  processes  to  a 
much  higher  level  of  fidelity  and  sophistication  than  in  the  Natick  effort.  Both  the 
STRICOM  and  Natick  efforts  would  have  to  have  some  minimal  level  of  performance  in 
both  the  structural  and  character  processes  in  case  both  efforts  are  not  funded.  That  is, 
both  efforts  need  to  be  stand-alone  products. 
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4.  CONCLUDING  REMARKS 

The  general  problem  addressed  by  the  first  phase  of  this  effort  was  developing  real¬ 
time  models  to  describe  weapon  effects  for  DIS  virtual  environments.  More  specifically, 
these  effects  needed  to  describe  building  structural  integrity  and  personnel  trauma  and 
be  suitable  for  implementation  in  a  DIS  virtual  environment.  The  original  Phase  I 
proposal  was  oriented  toward  large  weapons,  i.e.,  weapons  that  could  completely 
destroy  buildings  and  was  decoupled  from  the  virtual  representation  of  these  effects. 
That  is,  originally,  the  Phase  I  effort  was  conceived  of  as  theoretically  developing  the 
appropriate  weapon  effects  algorithms  which  would  then  be  implemented  in  Phase  II 
and  integrated  with  virtual  representations  outside  the  scope  of  this  effort. 

During  the  course  of  this  effort,  after  several  interchanges  with  USAIS,  it  became  clear 
that  the  most  relevant  effects,  at  least  initially,  were  small  arms.  This  change  of 
orientation,  in  a  sense,  was  fortunate  in  that  many  of  the  fundamental  problems 
associated  with  representing  weapon  effects  converge  in  modeling  small  arm 
interaction  with  building  materials  and  personnel.  Further,  we  decided  to  select  a 
subset  of  effects  (bullet  holes,  debris  and  ejecta,  and  residual  penetrators)  that  are 
common  to  all  weapon  effects  and  develop  models  for  these  effects  to  a  level  of 
maturity  beyond  what  was  originally  envisioned  in  the  Phase  I  effort.  This  was  done  in 
order  to  crystallize  the  really  technically  stressing  problems  associated  with  the  goal  of 
this  effort  and  to  highlight  what  really  needed  to  be  done  in  Phase  II. 

The  basic  problem  of  creating  a  hole  and/or  many  holes  in  walls,  along  with  debris  and 
ejecta,  and  describing  in  parallel  the  effect  of  residual  penetrators  exiting  the  wall  and 
striking  other  walls  and/or  people,  distills  many  of  the  technical  obstacles  that  need  to 
be  confronted  in  simulating  these  and  other  effects  using  physically  based  models  in  a 
virtual  environment  (as  opposed  to  scripts).  This  set  of  problems  is  fundamental  and  in 
a  sense,  the  ability  to  successfully  address  these  problems  represents  a  “litmus  test” 
indicating  whether  a  general  solution/approach  is  feasible. 

It  became  apparent  for  this  application  that  development  of  the  weapon  effects  models 
perse  was  not  technically  stressing.  Well  established  methodologies  and  data  bases 
exist  to  evaluate  blast,  penetration,  and  damage  to  buildings  and  personnel  to  varying 
levels  of  approximation.  At  MRC  we  have  a  plethora  of  these  models  and  data  bases 
some  of  which  we  have  been  funded  to  develop.  The  methodologies  selected  for  Phase 
II  implementation,  as  well  as  their  rationale,  are  discussed  in  the  Phase  II  proposal  for 
this  effort.  They  represent  well  established  and  validated  approaches.  In  particular,  the 
methodology  for  describing  weapon  effects  on  personnel  was  developed  by  MRC  and 
software  already  exists  to  address  this  problem. 

By  the  end  of  Phase  I  it  became  apparent  that  the  really  technically  stressing  issues 
concern  coupling/embedding  physically  based  models  in  the  virtual  environment  not 
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developing  real-time  algorithms  which  describe  these  effects.  This  orientation  is 
reflected  in  our  Phase  II  proposal. 

We  did  not  devote  significant  resources  in  Phase  I  in  areas  where  we  already  “had  the 
answer.”  It  is  recognized  however,  because  of  the  visualization  “choke”  points  that 
implementing  models  of  physical  processes  in  the  visual  environment  would  entail, 
trade-offs  needed  to  be  made  in  terms  of  the  detail  with  which  physical  processes  could 
be  modeled  and  the  degrees-of-freedom  present  in  interaction  with  the  virtual 
environment.  The  work  done  in  Phase  I  helps  to  identify  where  these  trade-offs  can  be 
made  without  affecting  the  utility  of  the  virtual  environment  as  a  training  and  mission 
rehearsal  tool. 


Mission  Research  Corporation 
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5.  Listing  of  MRSAW  Source  Code 


o  o  o  oo  ^  n  n  no  rfh  o  o  o  o  d^)  o  o  0^5  o  n 


$DEBUG 

C  -  set  up  interface  to  dos 

INTERFACE  TO  INTEGER* 2  FUNCTION  SYSTEM  [C]  (STRING) 

CHARACTER  STRING [ REFERENCE ] 

END 

^  ========================  end  interface  setup  ============================ 

C  PROGRAM  AIRBULL7 . FOR 

C  CREATED  AUGUST  1996  TO  DECSCRIBE  THE  MOTION  OF  VAROIUS  BULLETS  IN 

C  AIR.  BESIDES  THE  ACCELERATION  DUE  TO  GRAVITY,  THE  DECELERATION 

0  PRODUCED  BY  WIND  DRAG  IS  MODELLED  BY  A  FUNCTIONAL  FORM 

C  A*V**(2-N) 

WHERE  A  IS  THE  AIR  DRAG  COEFFICIENT  AND  N  IS  THE  AIR  DRAG  INDEX. 
THESE  VARIABLES  ARE  CONSTANT  WITHIN  REGIONS  OF  LOCAL  BULLET  VELOCITY 
WHICH  IS  ASSUMED  TO  STAY  IN  THE  DOMAIN  DURING  THE  FLIGHT  PATH.  THIS 
IS  VALID  FOR  TARGET  RANGES  DEPENDING  ON  MUZZLE  VELOCITY. 


C 

% 


INPUT  DESCRIPTION; 


VO  =  INITIAL  BULLET  SPEED  (MUZZLE  VELOCITY)  IN  FT/SEC 
DLTAD  =  BULLET  ORIENTATION  WITH  THE  HORIZONTAL  IN  DEGREES 
R  =  TARGET  RANGE  IN  FT.  =  NORMAL  DISTANCE  FOR  OBLIQUE  TARGET 
A  =  AIR  DRAG  COEFFICIENT 
N  =  AIR  DRAG  INDEX 

VWX  =  CROSSWIND  SPEED  IN  THE  HORIZONTAL  DIRECTION 
VWY  =  CROSSWIND  SPEED  IN  THE  VERTICAL  DIRECTION 
ALPA, BETA, GAMA  =  DIRECTION  COSINES  OF  THE  NORMAL  TO  THE 
TARGET  WITH  RESPECT  TO  THE  BULLET  XYZ  AXES  SYSTEM 
ER  =  COEFFICIENT  OF  RESTITUTION  FOR  BULLET-TARGET  IMPACT 
RIGG  =  MOMEMNT  OF  INERTIA  OF  THE  BULLET  ABOUT  ITS  CM 
XG  =  DISTANCE  OF  THE  CENTER  OF  MASS  FROM  THE  BULLET  TIP 
R1  =  DISTANCE  OF  FRONT  TIP  FROM  THE  CENTER  OF  MASS  G 
R2  =  DISTANCE  OF  REAR  TIP  FROM  THE  CENTER  OF  MASS  G 


ENERGY  LOST  DUE  TO  PENETRATION  :  ARRAY  ENRLST(4,4) 
GBL/MRC  EXPERIMENTAL  DATA  BASE 


ENRLST(I, J) , 


I  =  PROJECTILE  TYPE  ;  J 
=  1  FOR  SLUG 
=  2  FOR  lOOLEAD 
=  3  FOR  #04CU-PLA 
=  4  FOR  9MM 


TARGET  TYPE 

1  FOR  INTERIOR  WALL 

2  FOR  INTERIOR  DOOR 

3  FOR  FLOOR 

4  CEILING 


NUMBER  OF  PELLETS  (NP)  VS  RANGE  DISTRIBUTION  (R) 
NP=NPELL(NPRJ) *EXP(-ALFA* (R/RREF) * (BEETA-r/R) ) 


IMPLICIT  REAL*8 (A-H,0-Z) 

REAL* 4  RANO 

REAL*8  N,NN1,NN2,LB,LF 

REAL*8  L1,M1,N1,L2,M2,N2,L3,M3,N3,LP,MP,NP 
REAL*8  LXP,MXP,NXP,LZP,MZP,NZP 
REAL* 8  NC(7) 

INTEGER* 2  SYSTEM 
EXTERNAL  Y1D0T2 , Y2DOT2 


o  o  o  on 


CHARACTER* 15  TARGET ( 4 ) , PRO J ( 4 ) 

DIMENSION  ENRLST(4,4) ,NPELL(4) ,ALFA(4) ,BEETA(4) ,RF(4) 

DIMENSION  XT(500) ,YT(500) ,TP1X(500) ,TP1Y(500) ,TP1Z(500) , 
TP2X(500)  ,TP2Y(500) ,TP2Z(500)  , 

CMGX(500) ,CMGY(500) ,CMGZ(500) 

DIMENSION  VC(8),AC(7) 

COMMON/INITDAT/YIO , YIDO , Y20 , Y2D0 
COMMON/ FUNDAT/ A , N , VWX , VWY , PI 
COMMON/COMDAT/R 

COMMON/TRANS  DAT/ VIM , DLTAIM , XQ , YQ , ALPA , BETA , GAMA 
COMMON/SYSINFO/Ll , Ml , N1 , L2 , M2 , N2 , L3 , M3 , N3 

DATA  VC/3600.0,2600.0,1800.0,1370.0,1230.0,970.0,790.0,0.0/ 
DATA  NC/0. 45, 0.30, 0.0, -1.0, -3. 0,-1. 0,0.0/ 

DATA  AC/246.0,801.3 ,7599. 0,1. 045E7,1.578E13,1.685E7, 21386/ 

DATA  ENRLST/4.402E6,5.66722E6,5.24832E6,5.66722E6, 

6.41025E6,8.72053E6,11.76933E6,19.73025E6, 

6.9828E6,9.14848E6,14.20848E6,24.64473E6, 

0.4651E6,0.86822E6,0.5472188E6,0.70922E6/ 

DATA  TARGET/'  INT.  WALL','  INT.  DOOR','  FLOOR','  CEILING'/ 

DATA  PROJ/'  SLUG','  #00  LEAD', '#4  CU-PLA','9MM  BULLET'/ 

DATA  NPELL/1,16,27,1/ 

DATA  ALFA/0 . ODO , 116 . 774D0 , 13 1 . 94 6D0 , 0 . ODO/ 

DATA  BEETA/0 . ODO , 0 . 026D0 , 0 . 025D0 , 0 . ODO/ 

DATA  RF/0 . ODO ,120. ODO ,120. ODO ,0.0/ 

CONSTANTS 


PI=DATAN(1.0D0) *4. ODO 
DTR=PI/180.0DO 


DEFAULT  VALUES 
NDEBUG=0 


CALL  TLOCATE(7,0) 

CALL  CLS(2,0) 

I=SYSTEM( 'CLS'C) 

WRITE (*,*)  '  Executing  AIRBULL7  Code 

WRITE (*,*)  '  ' 


PAUSE  '  Please  enter  any  key  to  continue. . . . ' 


I=SYSTEM ( • TYPE  BULLSC. INF ' C) 

READ(*,*)  SC 
I=SYSTEM( 'CLS'C) 

WRITE (*,*)  '  Please  enter  bullet  weight  in  grams..' 
READ(*,*)  WTGM 

CONVERT  WEIGHT  IN  GRAINS 

WT=WTGM* 15.43236DO 

WRITE (*,*)  '  Please  enter  bullet  diameter  in  inches...* 
READ(*,*)  DI 


DETEREMINATION  OF  BALLISTIC  COEFFICIENT  BC  FOR  THE  ABOVE  BULLET 
THIS  RESULT  IS  EMPIRICAL:  FROM  MODERN  PRACTICAL  BALLISTICS:  PEJSA 

RNUM=(DI+0.50D0) /DI**2 
RDIN=415.0D0*(8-2.0D0*SC+SC**2) 

BC=WT*RNUM/RDIN 

CALL  TLOCATE(7,0) 

CALL  CLS(2,0) 


I=SYSTEM( 'CLS'C) 


WRITE(*,302)  SC 

FORMAT ('  Bullet  Shape  Class  =  ',f6.2,/) 

WRITE (*,304)  WTGM 

FORMATC  Bullet  Weight  ',G10.4,'  GM',/) 

WRITE (*,306)  DI 

FORMATC  Bullet  Diameter  ',G10.4,'  IN',/) 

WRITE (*,308)  BC 

FORMAT('  Computed  Ballistic  Coefficient  =  ',G10.4,/) 
OPEN (UNIT=4 , FILE= ' AIRBULL7 .IN') 

READ(4,*)  V0,DLTAD,R 


DETERMINE  MEYEVSKI  PARAMETERS  A  AND  N 
DO  322  1=1,7 

TST=(VC(I+1)-V0)/(VC(I)-V0) 

NCNT=I 

IF(TST.LE.O.ODO)  GOTO  324 
CONTINUE 

WRITE (*,*)  '  MUZZLE  VELOCITY  IS  OUTSIDE  DATA  RANGE..' 

WRITE (*,*)  '  PROGRAM  TERMINATED  ...' 

STOP 


324 


CONTINUE 


A=1.0D0/AC(NCNT) 

N=NC(NCNT) 

%  MODIFY  A  TO  INCLUDE  BALLISTIC  COEFFICIENT  BC  IN  INGALLS'S  TABLE 

A=A/BC 

C  READ (4,*)  A,N 

# 

WRITE(*,326)  A 

326  FORMAT('  Computed  Meyevski  Parameter  A  =  ',G10.4,/) 

WRITE (*,328)  N 

328  FORMAT ('  Computed  Meyevski  Parameter  N  =  ',G10.4,/) 

^  PAUSE  '  Please  enter  any  key  to  continue. . . . ' 

C  CALL  TLOCATE(7,0) 

C  CALL  CLS(2,0) 


I=SYSTEM( 'CLS'C) 

READ (4,*)  VWX,VWY 
C  READ (4,*)  ALPA, BETA, GAMA 

READ (4,*)  THZD,THYD 

h  COMPUTE  ALPA,  BETA  AND  GAMA 

THZR=THZD*DTR 

THYR=THYD*DTR 

.  ALPA=-DCOS (THZR) *DCOS (THYR) 

*  BETA=-DSIN(THZR) 

GAMA=-DCOS (THZR) *DSIN(THYR) 

READ(4,*)  ER 

READ(4,*)  RHO,BM,BR,LB,LF 
.  READ(4,*)  R1,R2 

'  READ(4,*)  RIGG,XG 


DLTAR=DLTAD*  DTR 

Y10=0.0D0 

Y20=0.0D0 

Y1DO=VO*DCOS (DLTAR) 
Y2D0=V0*DSIN (DLTAR) 


C 

c 


CROSSWIND  CORRECTION  OF  INITIAL  CONDITIONS 

Y1D0=Y1D0-VWX 

Y2D0=Y2D0-VWY 


C 


DLTAR=DATAN (YZDO/YIDO) 


on  ^  T)000 


DLTAD=DLTAR/ DTR 
C  END  CROSSWIND  CORRECTION 


TANDLTA=DTAN ( DLTAR) 
TMAX=20 . ODO*R/VO 
NT=500 


DT=TMAX/ (NT-1) 

OPEN (UNIT=8 , FILE= ' TYIEX . PLT ' ) 
OPEN (UNIT=9 , FILE= ' TY2EX . PLT ' ) 
OPEN (UNIT=10 , FILE= ' YIYIDEX . PLT ' ) 
OPEN (UNIT=11 , FILE= ' Y1Y2DEX . PLT ' ) 
OPEN (UNIT=12 , FILE= ' Y1Y2EX . PLT ' ) 

WRITE (8,*)  NT 
WRITE (9,*)  NT 
WRITE (10,*)  NT 
WRITE (11,*)  NT 


F0=V0**N/A 


NN1=N-1.0D0 
NN2=N-2 . ODO 

G  =  GRAVITATIONAL  ACCELERATION  IN  FT/S**2 
G=32.2D0 


I=SYSTEM( 'CLS’C) 

CALL  TL0CATE(7,0) 

CALL  CLS(2,0) 

WRITE (*,*)  '  Continue  Executing  AIRBULL7  Code  . ' 

WRITE (*,*)  '  ' 

WRITE(*,2)  R 
2  FORMATC  TARGET  RANGE  =  ',G10.3,'  FT  '  ,/) 

^  WRITE(*,4)  VO 

4  FORMATC  INITIAL  BULLET  SPEED  =  ',G10.3,'  FT/SEC  '  ,/) 

WRITE(*,6)  DLTAD 

6  FORMATC  INITIAL  BULLET  ORIENTATION  WITH  HORIZONTAL  ',F6.3,'  DEG 

^  C.  ',/) 

WRITE(*,47) 

47  FORMATC  CROSSWIND  SPEED:  ') 

WRITE (*,48)  VWX 

^48  FORMATC  HORIZONTAL  DIRECTION  =  ',G10.3,'  FT/SEC  ') 


r 


WRITE (*,52)  VWY 

52  FORMAT ('  VERTICAL  DIRECTION  =  ',010.3,'  FT/SEC  ') 
%  PAUSE  '  ENTER  ANY  KEY  TO  CONTINUE  . ' 


I=SYSTEM{ 'CLS'C) 

•C  CALL  TLOCATE(7,0) 

C  CALL  CLS(2,0) 

C  RESULTS  FROM  SHORT  DISTANCE  LINEAR  THEORY 

0  WRITE(*,8) 

8  FORMAT ('  RESULTS  FROM  SHORT  DISTANCE  LINEAR  THEORY...',/) 


C 

# 


c 

c 

c 


DO  10  1=1, NT 
T=(I-1)*DT 

Y1=(F0/N) *(l.OD0-(l.OD0-NNl*V0*T/F0) **(N/NN1) ) 

Y2=~(G*F0/(V0**2*NN2) ) * (Y1+ (FO/ ( 2 *NN1) ) * 

C  ( (1.0D0-N*Y1/F0) ** (2.0D0*NN1/N)-1.0D0) ) 

Y2  CORRECTION  FOR  SHOT  ANGLE  DLTA:  ADDED:  AKC 

Y2=Y2+Y1*TANDLTA 

Y1D=V0* ( 1 . 0D0-N*Y1/F0) ** ( 1 . ODO/N) 

Y2D=(Y1D*G*F0/(V0**2*NN2) ) * ( ( 1 . 0D0-N*Y1/F0) ** (NN2/N) -1 . ODO) 

VEL=DSQRT(Y1D**2+Y2D**2) 

CROSSWIND  CORRECTION  FOR  LOCATION 

Y1C0R=VWX*T 
Y2 COR=VWY *T- Y 1 * VWX/ VO 

Y1=Y1+Y1C0R 

Y2=Y2+Y2COR 

TH=ATAN{Y2D/Y1D) 

THD=TH/DTR 

TEST  TARGET  HIT  BY  THE  CONDITION  THAT 
ALPA*Y1+BETA*Y2+R=0 
HIT=ALPA*Y1+BETA*Y2+R 
IF(Yl.GE.R)  GOTO  14 
IF(HIT.LE.O.ODO)  GOTO  14 


WRITE(8,*)  T,Y1 
WRITE (9,*)  T,Y2 
WRITE (10,*)  Y1,Y1D 
WRITE(11,*)  Y1,Y2D 
WRITE (12,*)  Y1,Y2 

TSV=T 

Y1SV=Y1 

Y2SV=Y2 

VELSV=VEL 

THDSV=THD 

CONTINUE 


WRITE (*,12) 

FORMAT ( •  TARGET  NOT  YET  REACHED  BY  THE  BULLET. ...',/) 

GOTO  18 

CONTINUE 


WRITE(*,32) 

FORMAT ('  TARGET  HIT  LOCATION:  ' ,/) 

WRITE (*,34)  YISV 

FORMAT ('  HORIZONTAL  DISTANCE  =  ',G10.4,'  FT',/) 
WRITE (*,36)  Y2SV 

FORMAT ( ‘  VERTICAL  DISTANCE  =  ' , GIO . 4 , '  FT ' , / ) 


WRITE (*,38)  TSV 

FORMATC  TIME  TO  HIT  THE  TARGET  =  •,G10.4,'  SEC',/) 


WRITE (*,42)  VELSV 

FORMATC  IMPACT  VELOCITY  AT  TARGET  =  ',G10.4,'  FT/SEC ' 


WRITE (*,44)  THDSV 

FORMATC  BULLET  ORIENTATION  WITH  HORIZONTAL  =  ',G10.4, 


PAUSE  '  ENTER  ANY  KEY  TO  CONTINUE 


I=SYSTEM('CLS'C) 

CALL  TL0CATE(7,0) 
CALL  CLS(2,0) 


/) 

DEG' 


WRITE(*,22) 

FORMATC  RESULTS  FROM  EXACT  NONLINEAR  ANALYSIS  '  ,/) 


CALL  RKNLR2 ( Y1D0T2 , Y2D0T2 , DT, 0 . ODO , TMAX, IRR) 


PAUSE  '  Enter  any  key  to  continue . ' 

%  I=S YSTEM ( ' CLS ' C ) 

I=SYSTEM ( ' TYPE  PROJ . INF ' C) 

READ(*,*)  NPRJ 

01  GENERATING  3-D  PLOT  FILES  FOR  THE  BULLET  LOCATIONS 

OPEN (UNIT=7 , FILE= ' Y1Y2 . PLT ' ) 

OPEN (UNIT=5 , FILE= ' YIYID. PLT ' ) 

OPEN (UNIT=6,FILE='Y1Y2D. PLT ' ) 

^  NDAT=0 

DO  100  1=1,500 

READ(7,*,END=102)  XT(I),YT(I) 

READ (5,*)  DUM,RY1D 
READ (6,*)  DUM,RY2D 

^  NDAT=NDAT+1 

RVEL=DSQRT(RY1D**2+RY2D**2) 

CMGX(I)=XT(I) 

^  CMGY(I)=YT(I) 

CMGZ(I)=O.ODO 

C  FOR  BULLETS  ONLY 

IF (NPRJ. NE. 4)  GOTO  100 

*  TP1X(I)=XT(I)+R1*RY1D/RVEL 

TPIY ( I ) =YT (I ) +R1*RY2D/RVEL 
TP1Z(I)=0.0D0 


^  TP2X(I)=XT(I)-R2*RY1D/RVEL 

^  TP2Y(I)=YT(I)-R2*RY2D/RVEL 

TP2Z(I)=0.0D0 

100  CONTINUE 

^102  CONTINUE 

WRITE (*,202)  NDAT 

202  FORMAT ('  DATA  POINTS  FOR  PRE-IMPACT  MODE ',14,/) 

C  SAVE  PLOT  FILE  AND  TERMINATE  IF  NO  HIT  OCCURS . . . 

^  IF(VIM.EQ.O.ODO)  GOTO  674 

C  USER  PROMT  FOR  PENETRATION  OR  RICOCHET  MODE 

C  CURRENTLY  DATA  IS  NOT  AVAILABLE  FOR  DETERMINING  THIS  TRANSITION 

^  AUTOMATICALLY  BY  THE  CODE.  MORE  TESTS  ARE  NEEDED  TO  ESTABLISH 


THESE  THRESHOLD  SPEEDS  FOR  VARIOUS  TARGETS  AND  PROJECTILES 

I=SYSTEM( 'CLS'C) 

I=SYSTEM ( ' TYPE  TARGET . INF ' C) 

READ{*,*)  NTGT 
I=SYSTEM( 'CLS'C) 

WRITE (*,*)  '  YOU  HAVE  ENTERED  THE  PROJECTILE -TARGET  PAIR:' 

WRITE (*,404)  PROJ(NPRJ) 

FORMAT (  '  PROJECTILE :  ' , 15A , / ) 

WRITE (*,406)  TARGET (NTGT) 

FORMAT(  '  TARGET:  ',15A,/) 

VIN=VELSV 
WRITE (*,412)  VIN 

FORMAT ('  IMPACT  VELOCITY  =  ',G12.3,'  FT/SEC',/) 

THRESHOLD  VELOCITY  VPEN  FOR  PENETRATION 
ENTERING  TARGET  PENETRATION  MODE 
CALCULATE  VPEN  ENERGY  LOSS  FROM  DATABASE 
ELOSS=ENRLST (NPRJ , NTGT ) 

WEIGHT  IS  USED  INSTEAD  OF  MASS  AS  ENERY  LOSS  DATA  ARRAY 
ENRLST  HAS  BEEN  MODIFIED  ACCORDINGLY 

EIN=0 . 5D0*WTGM*VIN**2 
VPEN=DSQRT (2 . ODO*ELOSS/WTGM) 

WRITE (*,662)  VPEN 

FORMAT ('  THRESHOLD  VELOCITY  FOR  PENETRATION:  ',G10.3,'  FT/SEC',/) 
WRITE (*,*)  '  PLEASE  ENTER  THE  THRESHOLD  VELOCITY  FOR  PENETRATION' 

READ(*,*)  VPEN 
I=SYSTEM( 'CLS'C) 

DETERMINE  WHICH  MODE  WILL  OCCUR:  RICOCHET  MODE  OR  PENETRATION  MODE 

IF (VIN. LT. VPEN)  GOTO  408 

WRITE (*,*)  '  TARGET  PENETRATION  OCCURRED...' 


EIN=0 . 5D0*WTGM*VIN**2 


423 

C 

C 


c 


502 

C 


•5I6 

c 

# 

512 


VEX=VIN*DSQRT { 1 . ODO-ELOSS/EIN) 

I=SYSTEM( 'CLS'C) 

WRITE (*,422)  VIN 

FORMAT (/,'  TARGET  ENTRY  VELOCITY  =  *012 . 3 , 'FT/S ' ,/) 

WRITE (*,423)  VEX 

FORMAT ( '  TARGET  EXIT  VELOCITY  =  ' G12 . 3 , • FT/S • , / ) 

ALSO  ENTER 

TARGET  DAMAGE  AREA:  POINT  FOR  SLUGS  AND  BULLETS  (NPRJ=1,4) 

AREA  FOR  BUCKSHOT  (NTGT=2,3) 

NTST= (NPRJ-1) * (NPRJ-4 ) 

IF(NTST.NE.O)  GOTO  502 

ONE  SHOT  DAMAGE  FOR  SLUG  AND  BULLETS 

WRITE (*,  504)  Y2SV 

FORMAT ('  TARGET  HIT  AREA  =  ',G10.2,'  ABOVE  HORIZONTAL',/) 

GOTO  506 

CONTINUE 

FOR  BUCKSHOT  MODE  ONLY 

RUN  RANDOM  NUMBER  GENERATOR  FOR  PAIRS  OF  NUMBER  OF  PELLETS 
NPEL=NPELL ( NPRJ ) 

WRITE (*,516)  NPEL 

FORMAT ('  NUMBER  OF  PELLETS  HITTING  THE  TARGET:  ',13,/) 
SPRANG  =  SPRAY  ANGLE  IN  RADIAN 
SPRANG=DATAN ( BEETA (NPRJ) ) 

SPRANGD=SPRANG/DTR 
WRITE (*,512)  SPRANGD 

FORMAT ('  PELLET  SPRAY  ANGLE  =  ' , GIO . 3 , ' (DEG) ' ,/) 


^C  SPRAD=SPRAY  RADIUS  ON  THE  THE  TARGET 

S  PRAD=R*  BEETA ( NPRJ ) 

WRITE (*,611)  SPRAD 

•eil  FORMAT ('  SPRAY  RADIUS  ON  THE  TARGET  =  ',G10.3,/) 
PAUSE  '  ENTER  TO  CONTINUE....' 


I=SYSTEM( 'CLS'C) 


WRITE (*,518) 

518  FORMAT (•  HIT  LOCATIONS:',/) 

WRITE (*,522) 

522  FORMAT(T3,'  RADIAL  DISTANCE ', T25 , 'ANGULAR  LOCATION ' ,/) 


C 


# 


®C 

C 


# 

510 

^524 

C 


READ  ISEED  FOR  FROM  FILE  FOR  RANDOM  NUMBER  GENERATION 
OPEN (UNIT=11 , FILE= ' RANDOM . SED ' ) 

READ (11,*)  ISEED 
CLOSE (UNIT=11) 

OPEN (UNIT=11 , FILE= ' RANDOM . SED ' ) 

PELLET  DISTRIBUTION  VS  RADIAL  DISTANCE 
STORED  IN  PELLET. DIS 

OPEN (UNIT=14 , FILE= ' PELLET . DIS ' ) 

SMR1=0 . ODO 

DO  510  I=1,NPEL 

RNPEL=NPEL*1 . ODO 

RI=I*l.OD0 

FOLLOWING  EQUATION  RESULTS  FROM  THE  INVERSION  OF 
PELLET  DISTRIBUTION  FUNCTION  AS  A  FUNCTION  OF  R,r 

SMR2=R* (BEETA(NPRJ) -(RF(NPRJ)/R) *DLOG(RNPEL/RI) /ALFA(NPRJ) ) 

WRITE (14,*)  I,SMR2 

RAN=RAN0 (ISEED) 

RAD=SMR1+(SMR2-SMR1) *RAN 
ANG=360.0D0*RAN 

WRITE(*,524)  RAD,ANG 

SMR1=SMR2 

CONTINUE 

FORMAT (T5 , G12 . 4 , T22 , G12 . 4 ) 

REWRITE  ISEED  FOR  NEXT  USE 
WRITE (11,*)  ISEED 


•c  TARGET  TYPE 

C  PROJECTILE  TYPE 

C  TARGET  RANGE 

C  MUZZLE  VELOCITY 

^C  FINISH  TARGET  PENETRATION  MODE:  WRITE  PLOT  DATA  AND  EXIT 


on  o  ooo  o  n  no  oo 


GOTO  532 


408  CONTINUE 

r 

WRITE (*,*)  '  TARGET  PENETRATION  DID  NOT  OCCUR...' 

WRITE(*,204) 

204  FORMAT (/,'  ENTERING  PROJECTILE  REBOUND  MODE - ',/) 


PAUSE  '  ENTER  ANY  KEY  TO  CONTINUE 


I=SYSTEM( 'CLS'C) 

CALL  TLOCATE(7,0) 
CALL  CLS(2,0) 


ORIGIN  WILL  NOW  BE  SHIFTED  TO  Q  TO  DETERMINE  THE  PATH 
DUE  TO  BULLET  RICOCHETTING 


COORDINATES  CONVERTED  TO  USER  FRAME 

DIRECTION  COSINES  OF  NEW  AXES  SYSTEM  WITH  ORIGIN  AT  Q 

X2-AXIS:  (L1,M1,N1) 

Y2-AXIS;  (L2,M2,N2) 

Z2-AXIS:  (L3,M3,N3) 


CALCULATION  OF  IMPACT  KINETICS 

(ALPA, BETA, GAMA)  =  D.C.  OF  THE  NORMAL  TO  THE  TARGET  WRT. 
BULLET  X-Y  FRAME 

DLTAIMR=DLTAIM*  DTR 

Ll=DCOS (DLTAIMR) 

M1=DSIN ( DLTAIMR) 

N1=0.0D0 

WRITE (*,*)  '  DIRECTION  COSINES  OF  THE  IMPACT  DIRECTION;' 

WRITE(*,352)  L1,M1,N1 
^352  FORMAT (3G12. 5,/) 

WRITE (*,*)  '  DIRECTION  COSINES  OF  THE  IMPULSE  DIRECTION;' 

WRITE (*,352)  ALPA, BETA, GAMA 

C  INCIDENT  ANGLE  WITH  THE  IMPULSE  DIRECTION 

•  THI=DACOS (-L1*ALPA-M1*BETA-N1*GAMA) 

THID=THI/DTR 

WRITE (*,354)  VIM 

354  FORMATC  INCIDENT  VELOCITY;  ',  G12 . 5 ,' FT/SEC ',/ ) 


WRITE (*,356)  THID 

356  FORMAT (’  INCIDENT  ANGLE  ',010.5, 

C  '  DEG.  WITH  THE  IMPULSE  DIRECTION',/) 

WRITE (*,358)  ER 

^58  FORMAT('  Coeff.  of  Restitution  ',G12.5,/) 

TFAC=DTAN(THI) 

TH0=DATAN (TFAC/ER) 

TH0D=TH0/DTR 

WRITE (*,362)  THOD 

362  FORMAT (•  REBOUND  ANGLE  ',G10.5, 

C  '  DEG.  WITH  THE  IMPULSE  DIRECTION',/) 

VP=VIM*DSIN (THI) /DSIN (THO) 

^  WRITE (*,364)  VP 

364  FORMAT ('  REBOUND  VELOCITY:  ',G12.5,'  FT/SEC ' ) 


RIMP=BM*VIM*DSIN (THI+THO) /DSIN (THO) 


WRITE (*,366)  RIMP 

366  FORMATC  Impulse  =  '  ,G12. 5,  '  (GM. FT/SEC)  '  ,/) 
VRAT=VIM/VP 

.  FACIMP=RIMP/(BM*VP) 

'  LP=VRAT*L1+FACIMP*ALPA 

MP=VRAT*M1+FACIMP*BETA 
NP=VRAT*N1+FACIMP*GAMA 


CS=L1 *ALPA+M1*BETA+N1 *GAMA 
SN=DSQRT ( 1 . 0D0-CS**2 ) 

L3=(M1*GAMA-N1*BETA) /SN 
M3= (N1*ALPA-L1*GAMA) /SN 
N3=(L1*BETA-M1*ALPA)/SN 

L2=M3*NP-N3*MP 

M2=N3*LP-L3*NP 

N2=L3*MP-M3*LP 


C  DIRECTION  COSINES  OF  THE  NEW  AXES  SYSTEM  AT  Q 

WRITE (*,*)  '  DIRECTION  COSINES  OF  THE  NEW  AXES  SYSTEM  AT  THE  IMPAC 

^  CT  LOCATION ' 

WRITE(*,368)  LP,MP,NP 
368  F0RMAT(  /,'  LP,  MP,  NP  =  ',3G12.5) 


WRITE(*,372)  L2,M2,N2 
^372  FORMAT (  '  L2 ,  M2,  N2  =  ',3G12.5) 

WRITE(*,374)  L3,M3,N3 
374  FORMAT (  '  L3,  M3,  N3  =  ',3G12.5) 


TRANSFER  TO  BULLET  RICOCHET  MODE  OR  SPHERE  REBOUND  MODE 


o  n 


C  ACCORDING  AS  TH  VALUE  OF  NPRJ  EQUAL  TO  4  OR  1 

IF(NPRJ.EQ.4)  GOTO  544 
C  SPHERE  REBOUND  MODE 

C  ALSO  ENTER  SPHERE  REBOUND  MODE  CALCULATION 

C  RUN  RK2  FOR  DIFFERENT  MUZZLE  VELOCITY  AND  SHOT  ANGLE 

C  PERFORM  COORDINATE  TRANSFORMATION 

C  STORE  DATA 


I  THP=DACOS (MP) 

SNP=DSIN(THP) 

LZP=NP/SNP 

MZP=O.ODO 

NZP=-LP/SNP 

I 

LXP=-NZP 
MXP=0 . ODO 
NXP=LZP 

C  NEW  INITIAL  CONDITIONS 

Y10=0.0D0 
Y1D0=VP*SNP 
Y20=0.0D0 
Y2DO=VP*MP 

RESET  TARGET  DIRECTION  COSINES 
SET  TO  RETURN  TO  THE  SHOOTING  POINT 


ALPA=1.0D0 
BETA=0 . ODO 
GAMA=0 . ODO 


RM=R 


CLOSE (UNIT=7 ) 

CALL  RKNLR2 ( Y1DOT2 , Y2DOT2 , DT , 0 . ODO , TMAX , IRR) 

# 

C  GENERATING  3-D  PLOT  FILES  FOR  THE  SLUG  LOCATIONS 

OPEN (UNIT=7 , FILE= ' Y1Y2 . PLT ' ) 

OPEN (UNIT=5 , FILE= ' YIYID . PLT ' ) 

^  OPEN (UNIT=6,FILE='Y1Y2D. PLT ' ) 

NDATT=0 


DO  605  1=1,500 

READ(7,*,END=607)  XT(I),YT(I) 
READ (5,*)  DUM,RY1D 


READ (6,*)  DUM,RY2D 
NDATT=NDATT+1 

RVEL=DSQRT (RY1D**2+RY2D**2 ) 
^  CMGX(I)=XT(I) 

C  CMGY(I)=YT(I) 

C  CMGZ(I)=O.ODO 

605  CONTINUE 


607  WRITE (*,612)  NDATT 

612  FORMAT ('  DATA  POINTS  FOR  POST-IMPACT  MODE ',14,/) 


L1=LXP 

M1=MXP 

N1=NXP 


L2=0.0D0 

M2=0.0D0 

N2=0.0D0 

^  L3=LZP 

M3=MZP 
N3=NZP 

CALL  SYSCHG ( XT , YT , CMGX , CMGY , CMGZ , XGO , YGO , ZGO , NDATT , NDAT) 

%  BYPASS  BULLET  RICOCHET  MODE 

C  TRANSFER  TO  FINAL  FILE  WRITING  FOR  PLOT 

GOTO  532 

^44  CONTINUE 

C  BULLET  RICOCHET  MODE 

C  RESET  RM  TO  R  =  RANGE 

^  RM=R 

C  OPEN  AND  WRITE  DATA  FILE  FOR  THE  BULLET  RICOCHET  CODE  BULLPEN 


PAUSE  '  Enter  any  key  to  continue...' 

I=SYSTEM( 'CLS'C) 

C  CALL  CLS(7,0) 

C  CALL  TLOCATE(2,0) 

C  ENTERING  RICOCHET  MODE  AFTER  REBOUND 

C  RIGG  =  MOMENT  OF  INERTIA  ABOUT  CENTER  OF  MASS  G 

^  XG  =  DISTANCE  OF  C.G.  G  FROM  THE  BULLET  TIP 


(To  n  o  o  o  rro  a  ^  oo  no 


C  284  RIGG=RHO*(LB**3/12.0D0+LF*LB**2/4.0D0) 

284  CONTINUE 

THDN  =  BULLET  ANGULAR  VELOCITY  DUE  TO  IMPACT 
THDN=RIMP*XG*DSIN(THI) /RIGG 
OPEN {UNIT=2 , FILE= ' BULLRICO .IN') 


WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 
WRITE (2,*) 


RHO,BR 

LB,LF 

A,N,VWX,VWY 
•0.0' ,VP 
•0.0, 0.0' 
'90.0' ,THDN 
LP,MP,NP 
L2,M2,N2 
L3,M3,N3 
XQ,YQ,ZQ 
CONCl,CONC2 
'0.0' ,TMAX,DT 
NDEBUG 
RM 


2 


CONTINUE 


CLOSE (UNIT=2 ) 


WRITE (*,*)  '  Bullet  Ricochetting  Mode  after  Impact  ' 

I  =  SYSTEM ( 'BULLRICO 'C) 


WRITE (*,*)  '  Bullet  Ricochetting  Completed  ...' 

OPEN  DATA  FILE  CREATED  BY  BULRIC02.FOR  CODE 

I=SYSTEM( 'COPY  BXYl.PLT  BXY12.PLT'C) 

I=SYSTEM( 'COPY  BXGD.PLT  BXGD2.PLT'C) 

I=S YSTEM ( ' COPY  BXY2 . PLT  BXY2  2 . PLT ' C ) 

I=S YSTEM ( ' COPY  BXYG . PLT  BXYG2 . PLT ' C ) 

I=S YSTEM ( ' COPY  BXG . PLT  BXG2 . PLT ' C) 


OPEN (UNIT=9 , FILE= ' BXYl . PLT ' ) 
OPEN (UNIT=10 , FILE= ' BXY2 . PLT ' ) 
OPEN (UNIT=11 , FILE= ' BXYG. PLT ' ) 

READ(9,*)  NDUM 
READ (10,*)  NDUM 
READ (11,*)  NDUM 


ZQ=0.0D0 

XG0=XQ-L1*R1 

YG0=YQ-M1*R1 


ZG0=ZQ-N1*R1 


OPEN(UNIT=20 

, FILE= ' AIRBULL7 . DBG ' ) 

WRITE (20,*) 

•  D.C.  OF  NEW  COORDINATE  AXES  SYSTEM  AT  C.M 

WRITE (20,*) 

'  L1,M1,N1=' ,L1,M1,N1 

WRITE (20,*) 

'  L2,M2,N2=' ,L2,M2,N2 

WRITE (20,*) 

’  L3,M3,N3=' ,L3,M3,N3 

WRITE (20,*) 

'  C.M.  LOCATION  AT  IMPACT: ' 

WRITE (20,*) 

'  XG  =  ' , XGO 

WRITE (20,*) 

'  YG  =  ' , YGO 

WRITE (20,*) 

'  ZG  =  ',ZG0 

CHANGE  TO  USER  COORDINATE  SYSTEM  THROUGH  SYSCHG 


PAUSE  '  Enter  any  key  to  continue...’ 


I=SYSTEM( ’CLS'C) 

CALL  CLS(7,0) 
CALL  TLOCATE(2,0) 


READ  COORDINATES  OF  FRONT  TIP  IN  NEW  COORDINATE  SYSTEM 

NDATT=0 

NMAX=500 

DO  15  I=1,NMAX 
READ(9,*,END=16)  XT(I),YT(I) 

NDATT=NDATT+1 

CONTINUE 


WRITE (*,*)  '  Bullet  Ricochet  Phase  Data  Points  ’ ,NDATT 

0L1=L1 

0M1=M1 

0N1=N1 


L1=LP 

M1=MP 

N1=NP 

CALL  SYSCHG ( XT , YT , TPIX , TPl Y , TPIZ , XGO , YGO , ZGO , NDATT , NDAT ) 


READ  COORDINATES  OF  REAR  TIP  IN  NEW  COORDINATE  SYSTEM 


DO  25  1=1, NDATT 
READ(10,*)  XT(I),YT(I) 


*25  CONTINUE 

CALL  SYSCHG (XT , YT , TP2X , TP2 Y , TP2Z , XGO , YGO , ZGO , NDATT , NDAT) 

A 

C  READ  COORDINATES  OF  C.M.  IN  NEW  COORDINATE  SYSTEM 


DO  135  1=1, NDATT 
READ(11,*)  XT(I),YT(I) 
^35  CONTINUE 


CALL  SYSCHG ( XT , YT , CMGX , CMGY , CMGZ , XGO , YGO , ZGO , NDATT , NDAT) 


Ll=OLl 

Ml=OMl 

Nl=ONl 

NDAT=NDAT+NDATT 


CONTINUE 


^622 

62 


^45 

674 

C 


WRITE  FINAL  OUTPUT  FILE  OF  THE  LOCI  OF  TIP  1,  TIP2  AND  C.M. 
IN  THE  USER  COORDINATE  SYSTEM 

CONTINUE 

OPEN (UNIT=16 , FILE= ‘ CMG . PLT • ) 

IF(NPRJ.NE.4)  GOTO  622 

OPEN (UNIT=12,FILE='FTIP. PLT • ) 

OPEN (UNIT=14 , FILE= ' RTIP. PLT ' ) 

WRITE (12, 62)  NDAT 
WRITE (14, 62)  NDAT 

WRITE (16, 62)  NDAT 

FORMAT (14) 

DO  45  1=1, NDAT 

WRITE(16,*)  CMGX(I) ,CMGY(I) , CMGZ (I) 

IF(NPRJ.NE.4)  GOTO  45 

WRITE(12,*)  TPIX(I) ,TP1Y(I) ,TP1Z(I) 

WRITE(14,*)  TP2X(I) ,TP2Y(I) ,TP2Z(I) 

CONTINUE 

IF(NPRJ.NE.4)  GOTO  624 

WRITING  ALL  DATA  (TIP  1,  TIP  2,  G) 

OPEN (UNIT=18 , FILE= ' XYZALL. PLT ' ) 


J 


OPEN (UNIT=19 , FILE= ' TXYZ12 . PLT ' ) 

DO  50  I=1,NDAT 

WRITE(18,51)  TPIX(I) ,TP1Y(I) ,TP1Z(I) , 
C  TP2X(I) ,TP2Y(I) ,TP2Z(I) , 

C  CMGX(I) ,CMGY(I) ,CMGZ(I) 

CONTINUE 
FORMAT (9G12. 5) 


WRITING  TO  3D  GRAPHIC  FILE 
RANGE 

TARGET  NORMAL 

WRITE (19,*)  R 

WRITE ( 19 , * )  ALPA , BETA , GAMA 

WRITE (19, 62)  NDAT 


DO  55  1=1, NDAT 
T=I*DT 

WRITE (19, 51)  T,TP1X(I) ,TP1Y(I) ,TP1Z(I) , 
C  TP2X(I) ,TP2Y(I) ,TP2Z(I) 

CONTINUE 

GOTO  99 


CONTINUE 

OPEN (UNIT=18 , FILE= ' XYZCM. PLT ' ) 

OPEN (UNIT=19 , FILE= ' TXYZCM . PLT ' ) 

DO  630  1=1, NDAT 

WRITE (18, 51)  CMGX(I) ,CMGY(I) ,CMGZ(I) 
CONTINUE 

WRITING  TO  3D  GRAPHIC  FILE 
RANGE 

TARGET  NORMAL 

WRITE (19,*)  R 

WRITE ( 19 , * )  ALPA , BETA , GAMA 

WRITE (19, 62)  NDAT 

DO  635  1=1, NDAT 
T=I*DT 

WRITE (19, 51)  T,CMGX(I) ,CMGY(I) ,CMGZ(I) 


CONTINUE 


WRITE(*,24) 

FORMAT  ('  Exit  AIRBULL7  Code . ') 

STOP 

END 


SUBROUTINE  S YSCHG (XI , YI , XO , YO , ZO , XQ , YQ , ZQ , JN , NS ) 
IMPLICIT  REAL*8(A-H,L-Z) 

INTEGER* 4  NS 

DIMENSION  XI(500) ,YI(500) ,XO(500) ,YO(500) ,ZO(500) 


COMMON/SYSINFO/Ll , Ml , N1 , L2 , M2 , N2 , L3 , M3 , N3 
ZI=O.ODO 


KEEP  TOTAL  SUM  NS+JN  NOT  GREATER  ITS  ALLOCATED  LENGTH  OF  500 

JN1=500-NS 

NT=NS+JN 

IF (NT. GT. 500)  JN=JN1 
DO  10  1=1, JN 

X0(NS+I)=L1*XI(I)+L2*YI(I)+L3*ZI+XQ 
YO (NS+I ) =M1*XI ( I ) +M2  *YI ( I ) +M3  *ZI+YQ 
Z0(NS+I)=N1*XI(I)+N2*YI(I)+N3*ZI+ZQ 
CONTINUE 

RETURN 

END 


FUNCTION  Y1DOT2 (X,Y1,Y1D,Y2,Y2D) 
IMPLICIT  REAL*8 (A-H,0-Z) 

REAL* 8  N 

COMMON/FUNDAT/A , N , VWX , VWY , PI 
V=DSQRT(Y1D**2+Y2D**2) 


INDEX  N  FOR  RETARDATION  FUNCTION  SUCH  THAT  THE  DECELERATION 

IN  AIR  IS  GIVEN  BY 

DECL=A*V**(2-N) 


VREL  =  VELOCITY  OF  BULLET  RELATIVE  TO  AIR 


XREL=Y1D-VWX 

YREL=Y2D-VWY 

VREL=DSQRT (XREL**2+YREL**2 ) 

NO  CROSSWIND  FORM 
Y1D0T2=-A*V** ( 1 . ODO-N) *Y1D 

GAMA1=DATAN (XREL/YREL) 

GAMA=PI-GAMA1 

Y1D0T2=-A*VREL** ( 1 . ODO-N) *XREL 

RETURN 

END 

FUNCTION  Y2DOT2 (X,Y1,Y1D,Y2,Y2D) 
IMPLICIT  REAL* 8 ( A-H , O- Z ) 

REAL* 8  N 

COMMON/FUNDAT/A , N , VWX , VWY , PI 
V=DSQRT(Y1D**2+Y2D**2) 


INDEX  N  FOR  RETARDATION  FUNCTION  SUCH  THAT  THE  DECELERATION 

IN  AIR  IS  GIVEN  BY 

DECL=A*V**(2-N) 

G  =  ACCELERATION  DUE  TO  GRAVITY  IN  FT/S**2 
G=32.2D0 


VREL  =  VELOCITY  OF  BULLET  RELATIVE  TO  AIR 

XREL=Y1D-VWX 

YREL=Y2D-VWY 

VREL=DSQRT (XREL**2+YREL**2 ) 

NO  CROSSWIND  FORM 
Y2DOT2=-G-A*V** ( 1 . ODO-N) *Y2D 

GAMA=DATAN ( YREL/XREL) 

GAMA=PI-GAMA 

Y2DOT2=-G-A*VREL** ( 1 . ODO-N) *YREL 

RETURN 

END 


SUBROUTINE  RKNLR2 ( FUNl , FUN2 , HI , XI , XFF , IRR) 
IMPLICIT  REAL*8(A-H,K,0-Z) 


COMMON/INITDAT/YIO , YIDO , Y20 , Y2D0 
COMMON/COMDAT/R 


COMMON/TRANS DAT/ VIM , DLTAIM , XQ , YQ , ALFA , BETA , GAMA 


IRR=ERROR  CODE 


IRR=0 


DTR=DATAN ( 1 . ODO ) *4 . ODO/180 . ODD 

TEST  INPUT  DATA 

N={XFF-XI)/HI 
HI=(XFF-XI)/N 
IF (HI)  12,22,32 


IRR=1 


WRITE(6,21) 

FORMAT (•  NEGATIVE  INCREMENT  IN  RK2 ;  PROGRAM  HALTED') 
RETURN 


IRR=2 

WRITE(6,31) 

FORMAT ('  ZERO  INCREMENT  IN  RK2 ;  PROGRAM  HALTED') 

RETURN 

CONTINUE 

OPEN  OUTPUT  FILES 

OPEN (UNIT=2 , FILE= ' TYl . PLT ' ) 

OPEN (UNIT=3 , FILE= ' TY2 . PLT ' ) 

OPEN (UNIT=5 , FILE= ' YIYID . PLT ' ) 

OPEN (UNIT=6,FILE='Y1Y2D. PLT ' ) 

OPEN (UNIT=7 , FILE= ' Y1Y2 . PLT ' ) 

N1=N+1 

WRITE (2,*)  N 
WRITE (3,*)  N 
WRITE (5,*)  N 
WRITE(6,*)  N 
WRITE (7,*)  N 

INITIAL  VALUES:  TRANSFERRED  THROUGH  THE  COMMON  BLOCK  COMDAT 


INTEGRATION  OF  TWO  COUPLED  NONLINEAR  EQUATION 

Y1D2=FUN1(X,Y1,Y1D,Y2,Y2D) 

Y2D2=FUN2 (X, Y1 , YID, Y2 , Y2D) 


BY  FOURTH  ORDER  RUNGE-KUTTA  ALGORITHM 


SET  INITIAL  VALUES 


Y1=Y10 

Y2=Y20 

Y1D=Y1D0 

Y2D=Y2D0 


WRITE(2,*)  XI, Y1 
WRITE (3,*)  XI, Y2 
WRITE (5,*)  Y1,Y1D 
WRITE (6,*)  Y2,Y2D 
WRITE (7,*)  Y1,Y2 

HI2=HI/2.0D0 

HI8=HI/8.0D0 

X=XI 


DO  40  1=1, N 

CALCULATION  OF  Kl,  K2 ,  K3,K4 

K11=HI*FUN1(X,Y1,Y1D,Y2,Y2D) 
K12=HI*FUN2 (X,Y1,Y1D,Y2,Y2D) 

WRITE (*,*)  '  K11,K12: • ,K11,K12 

X1=X+HI2 

Y11=Y1+HI2*Y1D+HI8*K11 

Y21=Y2+HI2*Y2D+HI8*K12 

YlDl=YlD+Kll/2 . ODO 
Y2Dl=Y2D+K12/2 . ODO 

K2 1=HI*FUN1 (XI , Yll , YlDl , Y2 , Y2  D) 
K22=HI*FUN2(X1,Y1,Y1D,Y21,Y2D1) 

YlD2=YlD+K21/2 . ODO 
Y2D2=Y2D+K22/2 . ODO 

K31=HI*FUN1 (XI, Yll, Y1D2 , Y2 , Y2D) 
K32=HI*FUN2 (XI, Yl, YID, Y21, Y2D2) 

X2=X+HI 

Y12=Y1+HI*Y1D+HI2*K31 

Y22=Y2+HI*Y2D+HI2*K32 

Y1D2=Y1D+K31 

Y2D2=Y2D+K32 

K4 1=HI*FUN1 (X2 , Y12 , Y1D2 , Y2 , Y2D) 
K42=HI*FUN2 (X2 , Yl , YID, Y22 , Y2D2 ) 

NEXT  STEP  VARIABLES 

Y1=Y1+HI* (Y1D+ (K11+K21+K31) /6 . ODO) 


Y2=Y2+HI* ( Y2D+ (K12+K22+K32 ) /6 . ODO) 

NEXT  STEP  DERIVATIVES  OF  THE  VARIABLES 


Y1D=Y1D+ (Kll+2 . 0D0*K21+2 . 0D0*K31+K41) /6 . ODO 
Y2D=Y2D+ (K12+2 . 0D0*K22+2 . 0D0*K32+K42 ) /6 . ODO 

VEL=DSQRT(Y1D**2+Y2D**2) 

TH=DATAN ( Y2  D/ YID ) 

THD=TH/DTR 

X=X+HI 

TEST  TARGET  HIT  BY  THE  CONDITION  THAT 
ALPA*Y1+BETA*Y2+R=0 
HIT=ALPA*Y1+BETA*Y2+R 

IF(Yl.GE.R)  GOTO  42 
IF(HIT.LE.O.ODO)  GOTO  42 
WRITE  TO  FILES 


WRITE (2,*) 
WRITE (3,*) 

X,Y1 

X,Y2 

WRITE (5,*) 
WRITE (6,*) 
WRITE (7,*) 

Y1,Y1D 

Y1,Y2D 

Y1,Y2 

XSV=X 

Y1SV=Y1 

Y2SV=Y2 

VELSV=VEL 

THDSV=THD 

CONTINUE 

WRITE(*,24) 

FORMAT ( '  TARGET  NOT  YET  REACHED  BY  THE  BULLET  . . . ' , 

VIM=0 . ODO 
DLTAIM=0 . ODO 

GOTO  46 

CONTINUE 

TRANSFER  DATA  TO  THE  MAIN  ROUTINE  THORUGH  TRANSDAT 


XQ=Y1SV 

YQ=Y2SV 

VIM=VELSV 

DLTAIM=THDSV 


WRITE (*,26) 

FORMAT ('  TARGET  HIT  LOCATION:  ' ,/) 

WRITE(*,48)  YISV 

FORMAT (•  HORIZONTAL  DISTANCE  =  ',G10.4,'  FT',/) 

WRITE (*,52)  Y2SV 

FORMAT ('  VERTICAL  DISTANCE  =  ',G10.4,'  FT',/) 

WRITE (*,54)  XSV 

FORMATC  TIME  TO  HIT  THE  TARGET  =  ',G10.4,'  SEC',/) 

WRITE (*,56)  VELSV 

FORMATC  IMPACT  VELOCITY  AT  THE  TARGET  =  ',G10.4,'  FT/SEC',/) 


WRITE (*,58)  THDSV 

FORMATC  BULLET  ORIENTATION  WITH  HORIZONTAL  =  ',G10.4,'  DEG',/) 

CLOSE (UNIT=2 ) 

CLOSE (UNIT=3 ) 

CLOSE (UNIT=5) 

CLOSE (UNIT=6) 

CLOSE (UNIT=7) 


RETURN 

END 


THIS  SUBROUTINE  GENERATES  RANDOM  NUMBER 
LYING  BETWEEN  0  AND  1 
FROM  GIVEN  SEED  IDUM 


FUNCTION  RANG (IDUM) 

INTEGER  IDUM, IA,IM, IQ, IR, MASK 
REAL  RANG, AM 

PARAMETER  (IA=168G7 , IM=2147483647 , AM=1 ./IM, IQ=127773 , IR=2836 , 
*MASK=123459876) 

INTEGER  K 

I DUM=IEOR ( I DUM , MASK ) 

K=IDUM/IQ 

IDUM=IA* (IDUM-K*IQ) -IR*K 
IF  (IDUM.LT.G)  IDUM=IDUM+IM 
RANG=AM*IDUM 
IDUM=IEOR ( IDUM, MASK) 

RETURN 

END 


